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1.0  INTRODUCTION 


The  Nuclear  Weapons  Environments  Model  (NWEM)  code  development  project 
was  a  six  month  effort  to  create  a  nuclear  preprocessor  for  the  Air  Force 
Advanced  Missile  Computer  Model  (AMM)  being  developed  by  Logicon  Company. 

The  development  of  NWEM  was  viewed  as  an  evolutionary  process  where 
a  relatively  complete,  although  simplified,  set  of  launch  and  aimpoint 
exclusion  contour  methodologies  will  be  developed.  As  these  are  utilized 
in  NWEM  and  the  dominant  environments  for  the  scenarios  of  interest 
become  apparent,  the  most  important  can  be  improved  in  accuracy  and  new 
effects  or  damage  mechanisms  introduced.  The  NWEM  code  has  been  brought 
to  a  state  of  development  where  this  can  be  accomplished. 

1.1  OVERVIEW  OF  USERS'  GUIDE 

The  NWEM  Computer  Code  User's  Guide,  Volume  II  of  the  three  volume 
documentation  of  the  NWEM  project,  contains  the  description  of  the  code  for 
the  user.  In  Section  2  of  this  volume  the  code  is  described.  The  code 
structure  is  developed  and  explained.  The  purpose  and  function  of  each 
subroutine  is  discussed  together  with  a  list  of  input  and  output  variables, 
corimon  blocks  and  calls.  The  major  code  areas  and  subroutines  are  func¬ 
tionally  flowcharted.  The  variables  contained  in  each  common  block  are 
defined.  An  input  guide  and  sample  problem  description  are  presented. 
Section  3  contains  sufficient  information  to  allow  the  user  to  maintain 
and  update  the  code. 
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2.0  CODE  DESCRIPTION 


The  objective  of  the  NWEM  code  is  to  generate  land  based  booster  and 
land  targeted  ballistic  RV  exclusion  contours  which  result  from  a  single 
nuclear  detonation  of  an  incoming  RV  and  corresponding  vulnerability/hard¬ 
ness  criteria  for  the  booster  or  RV.  For  boosters,  NWEM  generates  exclu¬ 
sion  contours  for  blast,  dust,  pebble/ejecta,  thermal  radiation,  neutron, 
x-ray  and  high  altitude  EMP.  For  reentry  vehicles,  NWEM  generates  exclu¬ 
sion  contours  for  blast,  dust,  pebble/ejecta,  neutrons  and  gammas.  The 
philosophy  chosen  to  perform  the  computation  of  the  exclusion  regions  was 
one  of  simple  but  reasonably  accurate  models  of  the  environment  combined 
with  vulnerability  criteria  to  determine  detailed  exclusion  contours. 

These  detailed  contours  are  simplified  and  combined  to  give  exclusion  con¬ 
tours  which  satisfy  AMM  code  requirements.  The  advantage  of  explicitly 
performing  environment  calculations  for  generation  of  the  exclusion  con¬ 
tours  is  flexibility.  Changes  in  the  vulnerability  criteria,  trajectories 
or  vehicle  parameters  are  directly  reflected  in  the  contours;  therefore, 
a  direct  calculation  allows  a  wide  range  of  system  parameter  variation. 
Additionally,  the  existing  models  of  the  environments  can  be  replaced  with 
more  sophisticated  models  if  necessary  and  new  contours  easily  generated. 

The  NWEM  code,  at  its  present  state  of  development,  is  a  valuable 
tool  for  the  analysis  of  reentry  vehicle  fratricide  and  booster  launch  ex¬ 
clusion  regions.  The  code  is  quick  running,  uses  low  computer  core  require¬ 
ments  and  is  operating  on  both  IBM  and  CDC  hardware.  Two  modes  of  operation 
are  available.  A  detailed  output  option  can  be  used  by  the  nuclear  effects 
analyst  to  evaluate  relative  importance  of  the  environments,  effects  of 
uncertainties  and  sensitivity  to  vulnerability  level.  The  normal  output 
mode  is  designed  for  compatibility  with  the  AMM  targeting  and  systems  eval¬ 
uation  code. 

The  completed,  first  version  of  the  NWEM  code,  is  designed  to  test 
environment  generation  models,  system  utilization  options  and  system  vul¬ 
nerability  levels  as  they  affect  fratricide  and  launch  exclusion.  It 
was  not  designed  to  reproduce  the  best  calculations  of  environment  using 
state-of-the-art  codes;as  this  was  not  consistent  with  the  computer  run 
time  and  core  storage  requirements  or  with  the  program  schedule.  NWEM 
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should  be  viewed  as  a  parametric  analysis  tool.  In  this  role  it  can  be 
extremely  valuable  in  developing  increased  understanding  of  the  relative 
importance  of  competing  environments  and  effects  and  to  prioritize  areas 
for  research  in  environment  generation  and  systems  vulnerabil i ty. 

Because  extensive  efforts  were  made  to  program  in  American  Standard 
FORTRAN  and  use  structured  top-down  programming  practices,  the  NWEM  code 
has  had  only  minor  difficulties  being  made  operational  on  the  AF/PAC  IBM 
3032  and  TRW's  370/158  as  well  as  the  TRW  CDC  Cyber  174/7418  system  on 
which  it  was  developed.  Additionally,  to  make  the  transition  from  one 
computer  system  to  another,  care  has  been  taken  to  keep  the  input  and 
output  to  a  minimum.  The  only  input/output  files  used  by  the  code  are 
standard  input  (5)  and  output  (6)  with  an  additional  output  file  (16) 
for  the  AMM  exclusion  contour  data. 


The  various  core  size  and  time  requirements 


TRW  CDC  Cyber  174/7418 
TRW  IBM  360/158 
AF  PAC  IBM  3032 

2.1  PROGRAM  STRUCTURE 


Core 

100,000  Octal  Words 
200,000  Decimal  Bytes 
200,000  Decimal  Bytes 


are  as  follows: 

Run  Time 

60-120  Seconds 
120-180  Seconds 
30-60  Seconds 


The  NWEM  computer  code  calculates  booster  and  RV  exclusion  region 
contours  for  use  with  the  AMM  program.  NWEM  is  a  stand-alone  preprocesso 
of  relatively  small  core  size  and  short  running  time.  It  is  programmed 
in  American  Standard  FORTRAN  and  uses  structured  top-down  programming. 

The  driving  parameters  or  input  to  NWEM  are  the  vehicle  kill  criteria, 
trajectory  (boosters)  or  vehicle  dimension  parameters  (RVs),  and  the  nu¬ 
clear  weapon  characteri sties  and  transmission  data  for  the  detonating  RV. 
The  results  of  an  NWEM  run  are  the  AW  exclusion  region  descriptions  on 
the  AMM  output  file  (16)  and  detailed  contour  data  to  various  degrees  (se 
lectable  by  the  user)  on  the  standard  output  file  (6). 

Since  certain  environments  are  modeled  or  used  by  both  bonsters  and 
RVs,  a  philosophy  was  chosen  that  each  environment  and  its  related  ex¬ 
clusion  contour  would  be  calculated  sequentially  and  a  flag  used  to 
define  whether  the  calculation  is  for  a  booster  or  an  RV.  The  advantage 
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of  treating  the  environments  and  the  associated  contours  in  this  manner 
is  that  it  allowed  for  rapid  development/assembly  of  environment  models 
and  contour  generation  routines  since  some  of  the  routines  were  able  to 
be  assembled,  derived  or  modified  from  existing  programs. 

In  one  sense,  NWEM  is  really  a  collection  of  environment  exclusion 
region  generators  assembled  together  by  a  common  input  and  by  algorithms 
which  manipulate  their  outputs  into  final  contours  in  the  required  AMM 
format  and  within  the  AMM  restrictions.  Indeed,  this  philosophy  is  even 
carried  out  to  the  extent  that  the  radiation  environments  are  separated 
in  a  similar  form  from  the  mechanical  environments. 

Briefly,  the  functional  flow  of  NWEM  can  be  described  as:  (1)  ini¬ 
tialize  data  storage,  (2)  input  and  prepare  data  storage,  (3)  combine 
exclusion  contours  to  meet  AMM  requirements  for  radiation  and  mechanical 
environments,  (4)  construct  final  contours,  and  (5)  output  final  AMM 
exclusion  region  data.  For  an  MRV  problem,  an  additional  cal cul ational 
step  and  combination  of  multiple  burst  contours  occurs. 


2.2  SUBROUTINE  DESCRIPTION 


This  version  of  the  NWEM  computer  code  has  104  subroutines  or 
functions.  These  routines  vary  from  very  complex  to  simple.  To  provide 
the  user  with  a  basis  of  understanding  the  function  of  the  subroutines, 
a  format  for  presentation  has  been  selected  which  allows  the  routines 
to  be  decoupled  from  the  calling  routines  and  treated  individually.  For 
each  subroutine,  a  description  of  what  it  does,  the  argument  list  vari¬ 
ables  and  their  definitions,  the  variables  of  concern  from  each  common 
block  (defined  in  the  common  block  description),  data  statement  variables 
and  their  definitions  and  the  names  of  routines  called  are  given  if  ap¬ 
propriate.  Additionally  for  those  routines  which  control  the  cal culational 
flow  of  the  program,  functional  flow  diagram  are  included.  These  subrou¬ 
tine  descriptions  are  presented  on  the  following  pages  in  alphabetical 
order  except  for  the  main  routine  which  is  first. 


NWEM  Nuclear  Weapons  Environments  Model 

This  is  the  executive  routine  for  the  Nuclear  Weapons  Environments 
Model  computer  code  It  directs  the  order  of  input,  calculation 
and  output. 


Tape  file  5 


This  is  the  input  file 


Tape  file  6 
Tape  file  16 


This  is  the  program  output  file 

This  is  the  AMM  exclusion  region  data  output 
file 


Common  blocks  used 


/HTRNS/,  / ITRNS/ ,  /TRNSL/ ,  /TRNS/,  /XTITL/ 


IPASS 


Flag  indicating  number  of  passes  through  the 
calls  for  MRV's  of  the  third  type 


Routines  called: 


INTAL,  INPXT,  INPREP,  PTMASS,  ERR0UT,  0UTIN 
EMPN,  HYDR0,  RADIAT,  STATIC,  STATRV,  MRVSUM, 
AMM0UT,  AB0UT 


SUBR0UTINE  AB0UT  (NF) 


This  routine  stops  program  execution  either  for  a  normal  termination 
or  a  termination  resulting  from  code  error  messages. 

NF  input  Flat  to  indicate  type  of  termination 

(0/1  -  normal  /error) 

Routines  called:  ERR0UT 


SUBR0UTINE  AIR  (EEE,  RRR,  GM) 

This  routine  is  a  real  air  equation  of  state  model  based  on  work 
at  Air  Force  Weapons  Laboratory. 


EEE 

input 

Specific  internal  energy 

(ergs/g) 

RRR 

input 

Density  (g/cm3) 

GM 

output 

y-1  where  y  is  the  ratio 

of  specific  heats 

FUNCTI0N  AIRMAS(T) 

This  routine  defines  the  total  mass  loading  (kT)  of  the  dust  cloud. 


T 

input 

Time  after  burst  (sec) 

C0MM0N  /THREAT/ 

input 

I 

C0MM0N  /MASSY/ 

input 

Al,  A2 ,  A5 ,  A6,  Bl,  B2,  B5,  B6,  TO,  Tl,  T3, 
T4 ,  T5 ,  T6 ,  WD3,  WD4 ,  WD5 
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SUBROUTINE  ALTSET  (HI,  ALT,  NSAVE) 


This  routine  sets  up  the  altitude  break  points  used  for 
calculated  RV  trajectory. 

the 

HI 

input 

Initial  incoming  altitude  of  the 
{Maximum  used  is  300000.  ft) 

RV  (ft) 

ALT 

output 

Altitude  break  point  array  (ft) 

NSAVE 

input 

Numbe/  of  break  points  allowed. 

NSAVE -1 

will  be  calculated 


SUBROUTINE  AMM0UT 

This  routine  writes  the  AMM  exclusion  region  data  output  file  (16) 
and  if  selected  the  same  data  formatted  to  the  code  output  file  (6). 


C0MM0N  /HTRNS/ 

input 

All  variables 

C0MM0N  / ITRNS/ 

input 

ITYPE 

C0MM0N  /TRNS/ 

input 

RVANG 

C0MM0N  / R VC/ 

i  nput 

FR,  FRA 

C0MM0N  / BST/ 

input 

FBR,  FRC,  FRD 

C0MM0N  /XTITL/ 

input 

LEV0UT 

Routines  called: 

ERR0UT 

* 
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SUBR0UTINE  BF00T  (Tl,  ZO.  Z1 ,  Z2,  VDT,  YSTAR,  TSTAR,  IV,  DL,  WID,  CEN, 
RMAX,  RMIN ,  PHU,  PHI2,  IF00T) 


This  routine  constructs  an  entire  footprint  at  the  given  time  after 


burst  and 

then  defines 

the  simplified  footprint  dimensions. 

T 

input 

Time  after  burst  (sec) 

ZO 

input 

Ground  plane  altitude  (kft) 

Z1 

input 

Lead  RV  burst  altitude  (kft) 

12 

input 

Trailing  RV  burst  height  (RV)  or  missile 
launch  altitude  (booster)  (kft) 

YDT 

input 

Effective  burst  yield  for  blast  (Mt) 

YSTAR 

input 

Trajectory  offset  for  burst  above  sea 
level  (RV  only)  (kft) 

TSTAR 

input 

Trajectory  time  offset  for  burst  above 
sea  level  (RV  only)  (kft) 

IV 

input 

Vehicle  type.  2  for  booster,  4  for  RV 

DL 

output 

Rectangular  region  half-length  (kft) 

WID 

output 

Rectangular  region  half-width  (kft) 

CEN 

output 

Uprange/downrange  position  of  rectangle 
center  (kft) 

RMAX 

output 

Annular  region,  maximum  radius  (kft) 

RMIN 

output 

Annular  region,  minimum  radius  (kft) 

PHI1 

output 

Central  angle  of  uprange  (RV)  or  downrange 
(booster)  annular  segment  (degrees) 

PH  12 

output 

Central  angle  of  downrange  (RV)  or  uprange 
(booster)  annular  segment  (degree) 

IF00T 

output 

Flag.  IF00T  =  0  for  no  footprint 

C0MM0N  / XTITL/ 

input 

LEV0UT 

NMAX 

data 

Maximum  number  of  iterations  allowed  in 
determining  the  blast  footprint  (200) 

Y0LD 

data 

Internal  parameter  initialized  to  zero 

C0LO 

data 

Internal 

parameter 

initial ized 

to 

zero 

R0LO 

data 

Internal 

parameter 

initial ized 

to 

zero 

SIGN2 

data 

Internal 

parameter 

initial ized 

to 

zero 

SS 

data 

Internal 

parameter 

initial ized 

to 

zero 

SSAVS 

data 

Internal 

parameter 

initialized 

to 

zero 

SSN 

data 

Internal 

parameter 

initial ized 

to 

zero 

SSN2 

data 

Internal 

parameter 

initial ized 

to 

zero 

SSV 

data 

Internal 

parameter 

initial ized 

to 

zero 

Routines  called: 

FLY0UT, 

REENTR,  ERR0UT 

SUBROUTINE  BF00T  FLOW  DIAGRAM  (Con't) 


2  STEP 
SEARCH 


L 


YES- 


INTER°OLATE  FOR 
CRITICAL  POINT 


STORE  INTERCEPT 
POSITION 


/ FOOTPRINT 
\COMPLETE? 


SETUP  FOOTPRINT  FOLLOWER 


INCREMENT  X,  Y 

muz 


CALL  FLYOUT  OR  REENTER 


BLAST  INTERCEPT  CONDITIONS 


CRITICAL  INTERCEPT 


CONDITION  BRACKETED? 


h*-\ 


NO¬ 


CHANGE  SEARCH 
DIRECTION 


SUBROUTINE  BISRCH  (TMAX,  TMIN,  A,  TX,  TY,  XP,  YP,  ZP,  R,  RS,  XLP,  YLP,  T, 
YD,  Z1,  22,  IV) 

This  routine  performs  a  binary  search  to  determine  shock  front  and 
vehicle  intercept  time  and  position. 


TMAX 

input 

Maximum  limit  of  time  interval  to  be  searched 
(sec) 

TMIN 

input 

Minimum  limit  of  time  interval  to  be  searched 
(sec) 

A 

input 

Search  flag.  A  =  +1  for  missile  flying  into 
shock.  A  =  -1  for  missile  flying  out  of 
shock  region. 

TX 

output 

Intercept  time  after  burst  (sec) 

TY 

output 

Intercept  time  after  missile  launch  (sec) 

XP 

output 

Intercept  position  component  perpendicular 
to  trajectory  plane  (kft) 

YP 

output 

Intercept  position  horizontal  component 
parallel  to  trajectory  plane  (kft) 

ZP 

output 

Intercept  position.  Altitude  above  ground 
plane  (kft) 

R 

output 

Slant  range  from  burst  point  to  vehicle  (kft) 

RS 

output 

Shock  radius  (kft) 

XLP 

input 

Missile  launch  point  relative  to  burst. 
Component  perpendicular  to  trajectory,  (kft) 

YLP 

input 

Missile  launch  point  relative  to  burst. 
Component  parallel  to  trajectory  (kft) 

T 

input 

Launch  time  after  burst  (sec) 

YD 

input 

Effective  blast  yield  of  burst  (Mt) 

Z1 

input 

Burst  elevation  (kft) 

11 

input 

Launch  point  elevation  (kft) 

IV 

input 

Vehicle  type.  2  for  booster,  4  for  RV 

Routines  called: 

VL0C,  RSHK 
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SUBR0UTINE  BLAST  ( YD1 ,  ZO,  Z1 ,  Z2,  CR1,  CR2,  CR3,  IV) 


This  routine  is  a  sub-executive  which  controls  the  development  of 
the  time  dependent  blast  exclusion  region. 


YD1 

input 

Burst  yield  (Mt) 

ZO 

input 

Ground  plane  elevation  (kft) 

Z1 

input 

Lead  RV  burst  altitude  (kft) 

22 

input 

Trailing  RV  burst  height  (RV)  or  missile 
launch  altitude  (booster)  (kft) 

CR1 

input 

Blast  kill  criteria  one: 

Booster  -  overpressure  (psi) 

RV  -  total  acceleration  (g) 

CR2 

input 

Blast  kill  criteria  two: 

Booster  -  Overpressure/ ambient  pressure 
ratio 

RV  -  axial  acceleration  (g) 

CR3 

input 

Blast  kill  criteria  three: 

Booster  -  q  •  a  (psf-deg) 

RV  -  normal  acceleration  (g) 

IV 

input 

Vehicle  type:  2  booster,  4  for  RV 

C0MM0N  /BOAT/ 

output 

T.  BL,  BW,  BC,  BRMX,  BRMN ,  PHI,  PH2 ,  NR,  N( 

C0MM0N  /HTRAJ/ 

i  nput 

TIM,  RNG,  11,  ND 

C0MM0N  /CONST/ 

output 

THRO,  PI,  Cl,  Rl,  T1 

C0MM0N  / XTITL/ 

input 

LEV0UT 

Routines  called: 


MATM62 ,  GCHECK,  ERR0UT,  TERPL,  BF00T 


SUBROUTINE  CDBCMP  (COB) 


This  routine  calculates  the  base  drag  coefficient 


CD8 

output 

Coefficient  of  base  drag 

C0MM0N  /BC/ 

input 

TCDEG,  TCRAD,  SINTC,  C0STC,  RN,  RB 

C0MM0N  /FPSC/ 

input 

XM8 

C0MM0N  / GPC/ 

input 

GAM8 

C0MM0N  /TRNSC/ 

input 

LAMINR 

C0FW0N  /PTMC/ 

input 

PI 

TABM8 

data 

Mach  numbers  table 

tabco 

data 

Reference  base  drag  coefficients  corresponding 
to  the  Mach  numbers 

I  BO 

data 

Number  of  Mach  numbers  in  table 

Routines  called: 

FNFME,  FNFRE,  TERPL 
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SUBROUTINE  COFCMP  (CDP) 


5 


| 

I 

f 

I' 

; 

i 


( 


This  routine  calculates  the  coefficient  of  pressure  drag  on  the 
body. 


CDP 

output 

Coefficient  of  pressure  drag 

C0MM0N 

/BC  / 

input 

TCDEG,  SINTC,  C0STC,  RN,  R8 

C0MM0N 

/ FSPC/ 

input 

XM8 

C0MM0N  / GPC/ 

input 

GAM8 

XMTAB 

data 

Mach  numbers 

TAB06 

TAB09 

TAB01 2 

TAB01 5 

TAB26 

TAB29 

TAB212 

TAB215 

TAB46 

TAB49 

TAB412 

TAB415 

data 

Subsonic-transonic  forebody  pressure  drag 
coefficients  for  various  bluntness  ratios 
and  cone  hal f  angles 

(For  a  give  TABXZ  the  data  would 
be  bluntness  ratio  X,  cone  half 
angle  Z;  i.e.,  TAB412  is  for  a 
bluntness  ratio  of  4  and  cone 
hal f  angle  of  120) . 

TCTAB 

data 

Cone  hal f  angles 

BRTAB 

data 

Bluntness  ratios 

118 

data 

Number  of  Mach  numbers  and  associated 
subsonic-transonic  forebody  pressure  drag 
coefficients 

14 

data 

Number  of  bluntness  ratios  and  associated 
subsonic-transonic  forebody  pressure 
drag  coefficients 

Routines  called: 

TERPL,  QUAD,  FNFP 
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SUBR0UT I NE  COFLAM  (CDFL) 


This  routine  calculates  the  coefficient  of  laminar  frictional  drag 
on  the  body. 


CDFL 

output 

Coefficient  of  laminar  frictional  drag 

C0MM0N  / BC/ 

input 

TCDEG,  C0STC,  TANTC,  RN,  RB 

C0MM0N  /FSPC/ 

input 

U8 ,  RH08 ,  XMU8,  XM8,  REFT8,  Q8 

C0MM0N  / GPC/ 

input 

PR 

C0MM0N  / LPNC/ 

input 

All  variables 

Routines  called: 

CDFXQ,  FNFFL,  FNBLL 

SUBR0UTINE  CDFMIX 

(CDFTT,  CDF) 

This  routine  calculates  the  skin  fraction  drag  coefficient  when 
flow  is  mixed  laminar  and  turbulent. 

CDFTT 

input 

Turbulent  skin  friction  drag  coefficient 

CDF 

output 

Skin  friction  drag  coefficient 

C0MM0N  /BC/ 

input 

RN,  RB 

C0MM0N  /TRNSC/ 

input 

RNRTR 

Routines  called: 

CDFTRB,  CD FLAM 

SUBR0UTINE  CDFTRB 

(CDFT) 

This  routine 
on  the  body. 

computes 

CDRT 

output 

C0MM0N  / BC/ 

input 

C0MM0N  /FPSC/ 

input 

C0MM0N  / GPC/ 

input 

C0MM0N  /LPNC/ 

input 

Routines  cal  led: 

SUBR0UTINE  CDFXQ  (XM8,  TW) 

This  routine 

cal cula' 

number. 

XM8 

input 

TW 

output 

XMTAB 

data 

TWTAB 

data 

16 

data 

Routines  called: 

the  coefficient  of  turbulent  frictional  drag 

Coefficient  of  turbulent  frictional  drag 

TCDEG,  C0STC,  TANTC,  RN,  RB,  XK 
U8,  RH08,  XMU8,  XM8 ,  REFT8,  Q8 
PR 

All  variables 

CDFXQ,  FNFFT,  FNRGH,  FNBLT 


the  RV  body  temperature  for  the  given  Mach 

Mach  number 
RV  wall  temperature 
tech  numbers 

Body  wall  temperatures  corresponding  to  the 
Mach  numbers  (°R) 

Number  of  Mach  number  values  in  table 
TERPl 
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SUBR0UTINE  C0MP  (Tl,  T2,  N,  M,  IV) 


This  routine  produces  composite  footprints  durinq  a  given  time 
interval  based  on  maximum  dimensions  of  blast,  ejecta  and  dust 
footprint  data. 


Tl 

input 

Beginning  of  time  interval  (sec) 

T2 

input 

Ending  of  time  interval  (sec) 

N 

input 

Number  of  times  at  which  composite 
are  required 

footprints 

M 

input 

Index  which  defines  the  beginning 
for  storage  of  composite  footprint 
common  blocks  BST  and  RVC 

location 
data  in 

IV 

input 

Vehicle  type;  2  for  booster,  4  for 

RV 

C0MM0N  /BDAT/ 

input 

TB,  BL,  BW,  BC,  NR 

C0MM0N  /DDAT/ 

input 

TD,  DL,  DW,  DC,  ND 

C0MM0N  /EDAT/ 

i  nput 

TE,  EL,  EW,  EC,  NE 

C0MM0N  / BST/ 

output 

F0TLWC,  NF0 

C0MM0N  / R VC/ 

output 

TMFLWC,  NTM 

Routines  called: 


TERPL,  C0NVRT 


LOOP  OVER  PREVIOUSLY 
SELECTED  TIMES 


SUBROUTINE  C0MPAR  (IPHASE,  UMAX,  J) 


Calculates  maximum  radiation  contour  extent  over  all  times  for  each 
radiation  environment. 


IPHASE  input  See  /RAD/ 

UMAX  input  Maximum  number  of  times,  T(I),  for  which  lethal 

volume  size  is  calculated  for  the  current 
radiation  environment 

4  -  Surface  burst 

5  -  Lethal  volume  intersects  surface. 

Grazing  trajectory  is  below  the  earth's 
surface  so  PT1 7  places  the  last  point 
at  the  surface 

6  -  Lethal  volume  intersects  surface. 

SUBROUTINE  PTMAX  calculates  the  next 
to  last  data  point  (five  5) 

7  -  Final  time  corresponds  to  a  lethal 

volume  point  directly  below  the  burst 
.  point 


J 

input 

Index  of  environment  type 

1  -  Neutron  fluence 

2  -  Prompt  gatmia  peak  dose  rate 

3  -  X-ray  Fluence 

4  -  Thermal  radiation  fluence 

C0MM0N  /RRR/ 

input 

TIMAX,  TIMIN,  JINDEX,  IHMAX,  MAXITM,  tL00PS 

C0MM0N  /RRR/ 

input 

All  variables 

C0MM0N  /  RR1  / 

input 

RH 

C0MM0N  / RR1  / 

output 

RUP,  RDN 

C0MM0N  / XTITL/ 

input 

LEV0UT 

Routines  called: 

None 
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SUBROUTINE  C0NVRT  (HL,  CEN,  DU,  DD,  NP) 

This  routine  converts  rectangle  half-length  and  position  of  center 


to  uprange  and  down 

range  extend. 

HL 

i  nput 

Half-lengths  of  dynamic  rectangular  footprin 

CEN 

input 

Positions  of  rectangle  center  downrange  of 
burst  ground  zero 

DU 

output 

Uprange  extents  of  dynamic  rectangular  foot¬ 
print 

DD 

output 

Downrange  extents  of  dynamic  rectangular 
footprint 

NP 

input 

Number  of  points  in  dynamic  rectangle  arrays 

(Maximum  of  50) 


SUBR0UTINE  CR0SS  (XA,  XB,  NA,  NB,  TA,  TB,  XCRS,  TORS) 


This  routine  finds  the  intersection  point  of  two  lines  which  are 
defined  by  linear  interpolation  of  data  points. 


XA 

input 

Ordinate  data  for  first  line 

XB 

input 

Ordinate  data  for  second  line 

NA 

input 

Number  of  data  points  for  first  line 
(50  Maximum) 

NB 

input 

Number  of  data  points  for  second  line 
(50  Maximum) 

TA 

input 

Abscissa  data  for  first  line 

TB 

input 

Abscissa  data  for  second  line 

XCRS 

output 

Ordinate  of  intersection  point 

TCRS 

output 

Abscissa  of  intersection  point 

Routines  called: 

TERPL 
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SUBR0UTINE  DRAGX  (JTJF,  CDT0T) 


This  routine  is  a  sub-executive  routine  which  directs  the  computation 
of  the  drag  coefficients  used  in  calculating  the  RV  trajectory. 


JTJF 

input 

Trajectory  flag  indicating  that  this  is 
coupled  trajectory  drag  coefficient  cal 
culation 

CDTjOT 

output 

Total  drag  coefficient 

C0MM0N 

/DR  I  DR/ 

input 

ISTPA2 

C0MM0N 

/  DR  1 0  R/ 

output 

ITJF,  ISTPAS 

C0MM0N 

/ FS PC/ 

input 

HT,  XM8 

C0MM0N 

/  GPC/ 

i  nput 

GAM8 

C0MM0N 

/TRNSC/ 

input 

LAMINR 

C0MM0N 

/ TRN  SC/ 

output 

LAM IN R 

Routines  cal  led: 

SP,  LP,  REDBDY,  TRPTID,  T0NSET,  TRPTTJ , 
CD FLAM,  CDFTRB,  CDFMIX,  CDBCMP,  CDPCMP 

SUBR0UTIME  ORG  (AM, 

ALPHA, 

WA,  CA,  CN) 

This  routine 

interpolates  on  RV  drag  coefficient  data 

AM 

input 

RV  Mach  number 

ALPHA 

input 

RV  angle  of  attack  (degrees) 

WA 

output 

RV  weight  to  area  ratio  (lbs/ft2) 

CA 

output 

Axial  drag  coefficient 

CN 

output 

Normal  drag  coefficient 

C0MM0N  /DRAG/ 

input 

CNBP,  CXBP,  ALFBP,  AMBP ,  WN,  NA,  NM 

Routines  called: 

TERPL2 
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SUBR0UTINE  OUST  (WB,  H0B1,  H0B2,  ECRITX,  TQUITX,  TQUIT) 

This  routine  is  a  sub-executive  which  defines  cloud  avoidance  and 
dust  erosion  exclusion  regions. 


WB 

input 

Burst  yield  (Mt) 

H0B1 

input 

First  RV  burst  elevation  (ft) 

H0B2 

input 

Second  RV  burst  elevation  (ft)  or  booster 
launch  elevation  (ft) 

ECRITX 

input 

Critical  level  of  intercepted  kinetic  energy 
(ft-lb/ft2) 

TQUITX 

input 

Dust  cloud  cutoff  time  (sec) 

TQUIT 

input 

Pebble  avoidance  cutoff  time  (sec) 

C0MM0N  /BURST/ 

output 

HT,  HM,  RM,  RB,  RH0 ,  HBRV 

C0MM0N  /DIMEN/ 

input 

NY,  NX,  YMAX ,  DY ,  DX 

C0MM0N  /HTRAJ/ 

input 

RR,  AA,  NN 

C0MM0N  /DDAT/ 

output 

T,  HI,  HW,  YC,  NF 

C0MM0N  /XTITL/ 

input 

LEV0UT 

Routines  called: 


NUCLD,  SIZER,  TERPL ,  SL0PE ,  ENERGY,  ERR0UT 


SU8R0UTINE  DYC0MB  (IV) 


This  routine  manages  the  combination  of  blast,  ejecta  and  dust 
footprint  data  into  one  dynamic  footprint. 


IV 

input 

Vehicle 

type;  2 

for  booster,  4  for  RV 

C0MM0N  /BOAT/ 

input 

TB,  BL, 

BW,  BC, 

NR 

C0MM0N  /DDAT/ 

input 

TD,  DL, 

DC,  ND 

C0MM0N  / EDAT/ 

input 

TE,  EL, 

EC,  NE 

Routines  cal  led: 


C0MP ,  XLIM,  C0NVRT ,  CR0SS 


SUBROUTINE  DYC0MB  FLOW  DIAGRAM 
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SUBROUTINE  DYNAM  (IV) 


This  routine  controls  the  filling  of  static  and  dynamic  footprint 
arrays  for  blast,  ejecta  and  dust  footprints. 


IV 

input 

Vehicle  type;  2  for  booster,  4  for  RV 

C0MM0N  /BOAT/ 

input 

TB,  BL,  BW,  BC,  BRMX,  BRMN ,  PHI,  PH2 ,  NR,  NC 

C0MM0N  /STA/ 

output 

TMN ,  TMX ,  RHMX ,  RUPMN ,  RDNMX ,  I  DATA 

C0MM0N  / BST/ 

output 

BLANG,  NBL 

Routines  called: 

XLIM,  C0NVRT,  DYC0MB 

SUBROUTINE  DYNAM  FLOW  DIAGRAM 


i 


TIMES  BEFORE  BURST 
STORE  MAX.  BLAST  REGIONS 
IN  STATIC  FOOTPRINT  ARRAY 

I 

TIMES  AFTER  BURST 
STORE  BLAST  REGION  DIMENSIONS 
IN  BLAST  ARRAY 

,  T— - 

CALL  DYC0MB 


MANAGE  COMBINATION  OF 
BLAST,  EJECTA  AND  DUST 
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SUBR0UT INE  EARLY  (T,  HT,  HM,  RM,  RB) 

This  routine  defines  the  early  time  dust  cloud  dimensions  for  times 
after  burst  of  one  minute  or  less. 


T 

input 

Time  after  burst  (sec) 

HT 

output 

Dust  cloud  top  height  (kft) 

HM 

output 

Dust  cloud  middle  height  (kft) 

RM 

output 

Dust  cloud  middle  radius  (kft) 

RB 

output 

Dust  cloud  base  radius  (kft) 

C0W0N  /THREAT/ 

input 

I,  SH0LD,  H0B 

C0MM0N  /STUFF/ 

input 

AFB ,  AHM,  AHT,  ARM,  BFB ,  BHM, 
TB0  ,  TFS ,  TFW ,  TL0 

output 

RFB 

THIRD 

data 

One  third 

Routines  called: 


RBASE 


SUBR0UTINE  EJCLD  (T,  W,  SH0B ,  HT,  RM,  RH0) 

This  routine  defines  the  ejecta  cloud  dimensions  and  density- 


T 

input 

Time  afterburst  (sec) 

W 

input 

Burst  yield  (Mt) 

SH0B 

input 

Scaled  height  of  burst  (ft/kt1^3) 

HT 

output 

Ejecta  cloud  top  height  (kft) 

RM 

output 

Ejecta  cloud  maximum  radius  (kft) 

RH0 

output 

Ejecta  cloud  density 

C0MM0N  /EC0NS/ 

Routines  called: 

input 

BBB,  HIN,  HPK,  RIN,  RPK,  THP,  THT,  TRP,  TRT 

EJC0NS 

SUBR0UTINE  EJC0NS  (W) 

This  routine  defines  yield  dependent  constants  for  EJCLD. 

W  input  Burst  yield  (Mt) 

C0MM0N  /EC0NS/  output  BBB,  HIN,  HPK,  RIN,  RPK,  THP ,  THT,  TRP ,  TRT 


SUBROUTINE  EJECTA  (WB,  H0B1,  H0B2,  CRITD,  CRITN) 


This  is  a  sub-executive  routine  which  controls  definition  of  the 
ejecta  cloud  exclusion  region. 


WB 

input 

Burst  yield  (Mt) 

H081 

i  nput 

First  R V  burst  elevation  (ft) 

H0B2 

input 

Second  RV  burst  elevation  or  booster  launch 
elevation  (ft) 

CRITD 

input 

Critical  ejecta  particle  size  (dummy) 

CRITN 

input 

Critical  number  of  hits  (dummy) 

C0MM0N  /HTRAJ/ 

input 

AA,  RR,  NN 

C0MM0N  /EDAT/ 

output 

T,  HI,  HW.  YC,  NF 

C0MM0N  / XT ITL/ 

input 

LEV0UT 

Routines  called: 


EJCLD,  TERPL,  ERR0UT 


SUBROUTINE  EJECTA  FLOW  DIAGRAM 


START 


INITIALIZE 

PARAMETERS 


CALL  EJCLD 

EJECTA  CLOUD 
SIZE  &  DENSITY 


/  DENSITY  =  0  \ 
OR  2ND  RV  BURST 
\  ABOVE  CLOUD?  / 


COMPUTE: 

CLOUD  SHADOW 
RECTANGULAR 
FOOTPRINT 


DIAGNOSTIC 
PRINTOUT 
k  DESIRED?  / 


PRINT  EJECTA 
SUMMARY 


RETURN 
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SUBROUTINE  EMPN  (H0B,  Y  EMPl ,  ICT0R) 


This  routine  calculates  the  booster  exclusion  contour  for  the  high 
altitude  EMP  environment. 


H0B 

input 

Burst  height  (ft) 

Y 

input 

Burst  yield  (kt) 

EMPl 

input 

EMP  Criteria  -  peak  electric  field  strength 

(V/m) 

ICT0R 

output 

Flag  to  specify  existence  of  exclusion  contour 
(0/1  -  no  contour/ contour) 

C0MM0N  /  XTITL/ 

i  nput 

LEV0UT 

C0MM0N  /TRAJ/ 

i  nput 

T,  A,  N 

C0MM0N  / BST/ 

output 

F0RADI ,  F0TIMS 

RE 

data 

Radius  of  the  earth  (km) 

ALTL 

data 

Minimum  burst  altitude  for  EMP  environment 
model  (km) 

ALTC 

data 

Maximum  altitude  of  the  EMP  environment  (km) 

ALT(5) 

data 

Height  of  burst  data  for  EMP  environment  (km) 

YG1 (9) 

data 

Gamma  ray  yield  data  for  EMP  environment  (kt) 

EF(9,5) 

data 

Electric  field  strength  for  EMP  environment 
model  (V/m) 

Routines  called: 

ERROUT,  TERPL,  TERPL2 

FUNCTI0N  ENERGY  (X,  Z,  EM) 


This  routine  defines  intercepted  kinetic  energy  along  a  given  path 
through  a  given  dust  cloud. 


X 

input 

Initial  Vehicle  cross  range  coordinate  (kft) 

z 

input 

Initial  Vehicle  altitude  (kft) 

EM 

input 

Trajectory  slope  at  Z 

C0MM0N  /BURST/ 

input 

HT,  HM,  RM,  RB,  RH0,  HBRV 

C0MM0N  /HTRAJ/ 

input 

ALT,  VEL,  NN 

C3 

data 

Kinetic  energy  unit  conversion  constant 

EPS 

data 

Round  off  epsilon  (l.OxlO9) 

Routines  called: 

TERPL 

SUBR0UTINE  ERR0UT  (FMT,  X,  N,  NF,  NT) 

This  routine  outputs  the  data  to  either  the  AMM  output  file  (16) 
and  the  normal  output  file  (6). 


FMT 

input 

Array  containing  format  information  for 
outputing  data 

X 

input 

Array  containing,  the  data  to  be  output 

N 

input 

Number  of  data  values  to  be  output.  If  zero 
only  the  format  is  gutput 

NF 

input 

Flag  to  indicate  whether  this  output  should 
be  followed  by  termination  of  program 

NT 

input 

File  (6  or  16)  to  which  data  wil^  be  output 

Routines  called: 

AB0UT 
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FUNCTI0N  FENV  (RX,  RANGE,  J) 


Returns  the  nuclear  radiation  environment  value  for  a  point  in  space 
corresponding  to  a  given  areal  density  and  range  from  burst  using 
mass  integral  scaling. 

input  Air  mass  integral  (gm/cm2) 


RANGE 


input  Range  from  burst  (km) 


input 


C0MM0N  /0UT 


i  n  pu  t 


Routines  called: 


Index  of  environment  type  whose  environment 
is  defined  in  the  units  indicated 

1  -  Neutron  fluence  (neutrons/cm2) 

2  -  Prompt  gamma  peak  dose  rate  (rad(Si )/sec) 

3  -  X-ray  fluence  (cal /cm2) 

4  -  Thermal  radiation  fluence  (cal/cm2) 

All  variables 

Natural  logarithm  of  the  build  up  factor  for 
to  exoatmospheric  fluence.  This  has  been 
found  by  semi-log  interpolation  of  In 
(build  up  data)  as  a  function  of  air  mass 
integral  using  TERPL 

ibid,  for  peak  prompt  gamma  dose  rate 

ibid,  for  X-ray  fluence 

TERPL 


SUBR0UTINE  FLY0UT  (T,  XLP,  YLP ,  Z2  ,  ZO ,  Z1  ,  YD,  S,  IV,  IM) 


This  routine  controls 
waves  and  defines  the 

search  for  booster  intercepts  with  shock 
resulting  environment  severity. 

T 

i  nput 

Time  of  launch  after  burst  (sec) 

XL  P 

input 

Launch  point  distance  from  burst,  coordinate 
perpendicular  to  missile  trajectory  (kft) 

YLP 

input 

Launch  point  distance  from  burst,  coordinate 
parallel  to  missile  trajectory  (kft) 

12 

input 

Missile  launch  elevation  (kft) 

ZO 

input 

Ground  plane  elevation  (kft) 

Z1 

input 

Burst  elevation  (kft) 

YD 

input 

Effective  blast  yield  (Mt) 

S 

output 

Environment  severity.  Maximum  ratio  of 
environment  level  to  kill  criteria 

IV 

input 

Vehicle  type:  2  for  booster,  4  for  RV 

IM 

output 

Critical  environment  number  associated  v:ith  S 

C0MM0N  /HTRA 0/ 

input 

TIM,  RNG,  ALT,  NTRAJ 

Routines  called: 

TERPL,  RSHK,  BISRCH,  L0AD,  VL0C 
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SUBROUTINE  FLYOUT  FLOW  DIAGRAM 


SUBROUTINE  FNBLL  (BRATI0 ,  XMACA8 ,  THETAC,  BLL) 


This  routine  calculates  the  skin  friction  drag  blowing  correlation 
function  for  laminar  flow. 


BRATI0 

input 

RV  bluntness  ratio 

XMACH8 

input 

Mach  number 

THETAC 

input 

RV  cone  hal f  angle  (deg) 

BLL 

output 

Skin  friction  drag  blowing  correlation  in 
function  for  laminar  flow 

C 

data 

Laminar  skin  friction  drag  blowing  cor¬ 
relation  constants  for  slender  RVs 

D 

data 

Laminar  skin  friction  drag  blunt  RVs 

Routines  called:  SOLV 


SUBROUTINE  FNBLT  (BRATI0,  XMACH8 ,  THETAC,  BLT) 

This  routine  calculates  the  skin  friction  drag  blowing  correlation 
function  for  turbulent  flow. 


BRAT 10 

input 

RV  bluntness  ratio 

XMACH8 

input 

Mach  number 

THETAC 

input 

RV  Cone  half  angle  (deg) 

BLT 

output 

Skin  friction  drag  blowing  correlation 
function  for  turbulent  flow 

C 

data 

Turbulent  skin  friction  drag  blowing  cor¬ 
relation  constants  for  slender  RVs 

D 

data 

Turbulent  skin  friction  drag  blowing  cor¬ 
relation  constants  for  blunt  RVs 

Routines  called: 

S0LV 
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SUBROUTINE  FNFFL  (BRATI0,  XMACH8,  THETAC,  FFL) 


This  routine  calculates  the  smooth,  no  blowing  skin  friction  drag 
correlation  friction  for  laminar  flow. 


8RATI0 

input 

RV  bluntness  ratio 

XMACH8 

input 

Mach  number 

THETAC 

input 

RV  cone  half  angle  (deg) 

FFL 

output 

Smooth,  no  blowing  skin  friction  drag 
correlation  function  for  laminar  flow 

C 

da  ta 

Laminar  smooth,  no  blowing  skin  friction 
drag  correlations  constants  for  slender  RVs 

D 

data 

Laminar  smooth,  no  blowing  skin  friction  drag 
correlation  constants  for  blunt  RVs 

Routines  called: 

S0LV 

SUBROUTINE  FNFFT  (BRATI0, 

XMACH8,  THETAC,  FFL) 

This  routine 
correl ation 

calculates  the  smooth,  no  blowing  skin  friction  drag 
friction  for  turbulent  flow. 

BRAT 10 

i  nput 

RV  bluntness  ratio 

XMACH8 

input 

Mach  number 

THETAC 

input 

RV  cone  half  angle  (deg) 

FFT 

output 

Smooth,  no  blowing  skin  friction  drag 
correlation  function  for  turbulent  flow 

C 

data 

Turbulent  smooth,  no  blowing  skin  friction 
drag  correlation  constants  for  slender  RVs 

D 

data 

Turbulent  smooth,  no  blowing  skin  friction 
drag  correlation  constants  for  blunt  RVs 

Routines  called: 

S0LV 
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SUBR0UTINE  FNFME  (BRATI0,  XMACH8,  THETAC,  FME) 


This  routine  calculates  the  local  Mach  number  correlation  function. 


BRAT10 

input 

RV  bluntness  ratio 

XMACH8 

input 

Mach  number 

THETAC 

input 

RV  cone  half  angle  (deg) 

FME 

output 

Local  Mach  number  correlation 

function 

C 

data 

Local  Mach  number  correlation 

constants 

Routines  called: 

S0LV 

SUBR0UTINE  FNFP  (BRATI0 ,  XMACH8,  THETAC,  FP) 

This  routine  calculates  the  forebody  pressure  drag  correlation 
function . 


BRAT 10 

input 

RV  bluntness  ratio 

XMACH8 

input 

Mach  number 

THETAC 

input 

RV  cone  half  angle  (deg) 

FP 

output 

Forebody  pressure  drag  correlation 

function 

C 

data 

Forebody  pressure  drag  correlation 

constants 

Routines  called: 

S0LV 

45 


SUBR0UTINE  FNFRE  (BRATI0,  XMACH8,  THETAC,  FRE) 

This  routine  calculates  the  local  Reynolds  number  correlation 
function . 


BRATI0 

input 

RV  bl  untness  ratio 

XMACH8 

input 

Mach  number 

THETAC 

input 

RV  cone  half  angle  (deg) 

FRE 

output 

Local  Reynolds  number  correlation 

function 

C 

data 

Local  Reynolds  number  correlation 

constants 

Routines  called: 

S0LV 

SUBR0UTINE  FNRGH  (BRATI0,  XMACH8,  THETAC,  RGH) 

This  routine  calculates  the  skin  friction  drag  roughness  correlation 
function . 


BRATI0 

input 

RV  bluntness  ratio 

XMACH8 

input 

Mach  number 

THETAC 

input 

RV  cone  hal f  angle 

(deg) 

RGH 

output 

Skin  friction  drag 

roughness  correlation 

C 

data 

Skin  friction  drag 
slender  RVs 

roughness  correlation 

for 

D 

data 

Skin  friction  drag 

roughness  correlation 

for 

blunt  RVs 

Routines  called:  S0LV 
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SUBR0UTINE  FSP 


This  routine  calculates  the  free  stream  properties  for  RV  trajectory 
drag  calculations. 


C0MM0N  /FSPC/ 

input 

U8,  P8, 

T8 ,  RH08,  XM8,  Q8 

C0MM0N  / FSPL/ 

output 

H8 ,  HT, 

XMU8, 

REFT8,  P8PT8 

C0MM0N  / GPC/ 

input 

CP8,  G, 

GAM8, 

XJ 

SUBR0UTINE  FSPX 


This  routine  calculates  ambient  free  stream  properties  for  the  RV 
trajectory  from  ambient  atmosphere  properties  and  RV  velocity. 


C0MM0N  / FSPC/ 
C0MM0N  /FSPC 

Routines  called: 


input  ALT,  U8 

output  P8,  T8»  A8,  RH08 ,  XM8,  Q8 

MATM62 


SUBR0UTINE  GCHECK  (GMAX) 

This  routine  finds  the  maximum  acceleration  along  normal  reentry 
trajectory. 


GMAX 

C0MM0N  /HTRAJ/ 
Routines  called: 


output  Maximum  acceleration  along  RV  trajectory  ( g ' s ] 
input  ALT,  VEL,  ND 
MATM62,  DRG 


i 


SU8R0UTINE  HPREP  (WK,  HI,  H2,  DCR1,  DCR2,  WM,  ZO,  Z1 ,  Z2,  ECRIT,  (QUIT) 

This  routine  zeros  footprint  data  arrays  and  converts  input  data 
units. 


WK 

input 

Burst  yield  (kt) 

HI 

input 

Lead  RV  burst  height  (ft) 

H2 

input 

Trailing  RV  burst  height  or  booster  launch 
site  elevation  (ft) 

0CR1 

input 

Dust  erosion  criteria  one  (kJ/cm2) 

0CR2 

input 

Dust  erosion  criteria  two  (min) 

WM 

output 

Same  as  WK  (Mt) 

ZO 

output 

Ground  plane  elevation  (kft) 

Z1 

output 

Same  as  HI  (kft) 

Z2 

output 

Same  as  H2  (kft) 

ECRIT 

output 

Same  as  DCR1  (ft-lb/ft2) 

TQU1T 

output 

Same  as  DCR2  (sec) 

C0MM0N  /TRAJ/ 

input 

Tl,  Rl,  Al,  VI,  AN1 ,  N1 

C0MM0N  /HTRAJ/ 

output 

T2 ,  R2 ,  A2 ,  V2,  AN2 ,  N2 

C0MM0N  /BOAT/ 

output 

TB,  BL,  BW,  BC,  BRMX,  BRMN,  PHI,  PH2 , 

NR,  NC 

C0MM0N  /DDAT/ 

output 

TO,  OL,  DW,  DC,  NO 

C0MM0N  /EDAT/ 

output 

TE,  EL,  EW,  EC,  NE 

C0MM0N  /STA/ 

output 

TMN,  TMX,  RAMX,  RUPMN ,  RDNMX,  IDATA 

SUBR0UTINE  HYDR0  (WK,  HI,  H2 ,  BCR1,  BCR2,  BCR3,  DCR1,  DCR2 ,  PCR1,  PCR2 , 
IPHASE) 

This  is  a  sub-executive  routine  which  drives  BLAST,  EJECTA,  DUST, 

and  DYNAM. 


WK 

input 

Burst  yield  (kt) 

HI 

input 

Lead  RV  burst  height  (ft) 

H2 

input 

Trailing  RV  burst  height  or  booster  launch 
site  elevation  (ft) 

BCR1 

input 

Blast  kill  criteria  one 

Booster  -  ratio  of  overpressure  to 
ambient  pressure 

RV  -  total  acceleration  (g) 

BCR2 

input 

Blast  kill  criteria  two 

Booster  -  dynamic  pressure  times  angle 
of  attack  (psf-deg) 

RV  -  axial  acceleration  (g) 

BCR3 

input 

Blast  kill  criteria  three 

Booster  -  dynamic  times  angle  of  attack 
(psi-deg) 

RV  -  normal  acceleration  (g) 

DCR1 

input 

Dust  erosion  criteria  one,  intercepted 
kinetic  energy  (kj/cm2) 

DCR2 

input 

Dust  erosion  criteria  two,  cloud  cut-off 
time  (min) 

PCR1 

i  n  pu  t 

Pebble/ejecta  penetration  criteria  one, 
critical  particle  diameter  (cm) 

PCR2 

input 

Pebble/e jecta  penetration  criteria  two, 
critical  hit  number  density  (#/cm2) 

IPHASE 

input 

Vehicle  type;  2  booster,  4  for  RV 

Routines  called: 


HPREP,  BLAST,  EJECTA,  DUST,  DYHAM 


FUNCTION  IBY  (UX,  X,  N) 

This  function  performs  a  binary  search  of  a  data  array  for  the 
index  of  lower  data  point  for  a  given  value. 

UX  input  Given  value  for  which  index  is  needed 


input  Data  array  to  be  searched 
input  Number  of  data  in  the  array 


SUBROUTINE  INLINE 


Merges  blast  exclusion  region  into  the  dust/ejecta  region  as  a 
time-dependent  set  of  rectangular  exclusion  contours.  Used  for 
MRV  options  2  or  3  during  the  boost  phase  only.  AMM  units  are 
used  throughout. 


C0MM0N  / BST/ 

input  BLANG,  F0TLWC,  NBL ,  NF0 

IIIMB0 

Flag  indicating  the  status  of  fitting  the 
blast  contour  times  TMINIT(I)  into  the  dy¬ 
namic  contour  time  array,  F0TLWC(1,J): 

0  -  Previous  time  has  been  evaluated 
select  the  next  blast  time  for 
merging  blast  data  into  the  dynamic 
array 

1  -  Continue  searching  for  the  dynamic 
contour  time  interval  in  which  the 
blast  contour  time  is  located 

TMINIT( 7) 

Blast  contour  times  with  units  changed  to 
conform  with  the  dynamic  times  (min) 

RBLAST (7) 

Blast  radius  at  times  corresponding  to  the 
dynamic  contour  times  (NMi) 

RCTDN( 7 ) 

Maximum  extent  of  the  final  exclusion  contour 
in  the  downrange  direction  relative  to  burst 
ground  zero  at  each  time  (NMi) 

RCTUP(7) 

Minimum  extent  of  the  final  exclusion  contour 
in  the  uprange  direction  relative  to  burst 
ground  zero  at  each  time  (NMi) 

WBMAX 

Maximum  radius  of  the  blast  contour  within 
the  current  dynamic  contour  time  interval  (NMi  ) 

TBMAX 

Time  corresponding  to  WBMAX  (min) 

Routines  called: 


TERPL 


SUBR0UTINE  INPREP  ( IPASS) 

This  routine  takes  the  input  data  and  transfers  the  data  to  the 
appropriate  working  storage  of  the  program. 


IPASS 

input 

Flag  indicating  pass  number  for  MRV  third 
type 

C0MM0N  / ITRNS/ 

input 

ITYPE,  NTRJ ,  NTRJX,  NDG1 ,  NDG2,  INEUT, 
INET,  IXRAY,  IXR 

C0MM0N  /TRNS/ 

input 

DUST!  ,  DUST2 ,  THERE,  XGM1 ,  YLD,  XNEUTX, 
XFRAC,  GAMDS,  GAMPL,  TIMX,  ALT,  RNG,  VEL, 
ANG,  AMX,  ANG,  ANGX,  CN,  CA,  RH0RN,  XNETX, 
RH0RX,  XPHI 

C0MM0N  /TRNS/ 

output 

DUST3,  THER1 ,  XTM1 ,  XTM2 

C0MM0N  /DRAG/ 

output 

XCN,  XCA,  XA1 ,  XAM,  NA,  NM 

C0MM0N  /0UT/ 

output 

All  variables 

C0MM0N  /TRAJ/ 

output 

All  variables 

C0MM0N  / WPN/ 

output 

All  variables 

C0MM0N  /S0R/ 

input 

All  variables 

C0MM0N  /TRAJX/ 

input 

All  variables 

C0MM0N  /THER/ 

input 

All  variables 

C0MM0N  /XTITL/ 

input 

LEV0UT 

LEVTYP(5) 

data 

Alphanumeric  representation  of  characters 
(blank,  0,  1,  2  and  3)  for  output  option 

LEVTYX(5) 

data 

Numeric  values  associated  with  LEVTYP 
(0,  0,  1 ,  2  and  3) 

Routines  called: 

ERR0UT,  TNFR 
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SUBROUTINE  INPXT  (IPASS) 


This  routine 

reads  the 

problem  input  to  the 

program. 

IPASS 

i  nput 

Flag  indicating  pass 
type 

number  for 

MRV  third 

C0MM0N  /HTRNS/ 

output 

All  variables 

C0MM0N  / ITRNS/ 

output 

All  variables 

C0MM0N  /TRNSL/ 

input 

All  variables 

C0MM0N  /TRNS/ 

output 

A1 1  variables  except 

DUST3,  XTM1 

and  XTM2 

C0MM0N  / XTITL/ 

output 

All  variables 

Routines  called: 

ERR0UT 

SUBROUTINE  INTAL 

This  routine  sets  the  values  of  the  input  common  block  storage 
and  AMM  output  common  block  storage  to  zero  (0)  to  start  the 


problem. 

C0MM0N  /HTRNS/ 

output 

All 

variables 

C0MM0N 

/ITRNS/ 

output 

All 

variables 

C0MM0N 

/TRNSL/ 

input 

All 

variabl es 

C0MM0N 

/TRNS/ 

output 

All 

vari  abl  es 

C0MM0N 

/  R  VC/ 

output 

All 

variables 

C0MM0N 

/  BST/ 

output 

All 

variables 
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SUBROUTINE  J0UT  (TM,  RAD,  RADI ,  RAD2 ,  J,  CRIT,  1PHASE) 


This  routine  outputs  the  radiation  exclusion  region  data  for  each 
radiation  environment. 


TM 

input 

Time  extent  of  the  radiation  exclusion 
region  (sec) 

RAD 

input 

Radial  extent  of  the  static  exclusion  region 
for  boosters.  Half  width  of  static  exclusion 
region  for  RVs  (km) 

RADI 

input 

Center  of  the  static  exclusion  region  for 

RVs  (km) 

RAD2 

input 

Half  length  of  static  exclusion  region  for 

RVs  (km) 

J 

input 

Index  of  radiation  environment 

CRIT 

input 

Radiation  environment  criteria  array 

l PHASE 

input 

Flag  to  indicate  RV  or  booster  exclusion 
region 

ICALL 

data 

Header  flag 

Routines  called: 

ERR0UT 

55 


p»*<V  «*•**<* 


■  /«ri* 


1 

SUBR0UTINE  LATE  (T,  HT,  HM,  RM,  RB) 

This  routine  defines  the  late  time  dust  cloud  dimensions  for  times  - 

after  burst  greater  than  one  minute.  j 


T 

input 

Time 

after 

burst  (sec) 

HT 

output 

Oust 

cloud 

top  height  (kft) 

HM 

output 

Dust 

cloud 

middle  height  (kft) 

RM 

output 

Oust 

cloud 

middle  radius  (kft) 

RB 

output 

Dust 

cloud 

base  radius  (kft) 

CpMMpN  /THREAT/ 

input 

I 

C0MM0N  /STUFF/ 

input 

ARM, 

BRM,  ! 

iMP,  HMS ,  HMO,  HTP , 

TPH, 

TRS 

Routines  called: 

RBASE 

SUBROUTINE  LOAD  (TB,  RS ,  YO,  X,  Y;  Z,  TSTAR ,  Z1  ,  IV,  S,  IM) 


This  routine  computes  blast  environments  at  shock  intercept  and 
returns  the  maximum  severity. 


TB 

input 

Time  before  impact  (RV)  or  after  launch 
(booster)  (sec) 

RS 

input 

Shock  front  radius  (kft) 

YD 

input 

Effective  blast  yield  (Mt ) 

X 

input 

Vehicle  x  coordinate  relative  to  burst  (kft) 

Y 

input 

Vehicle  y  coordinate  relative  to  burst  (kft) 

Z 

input 

Vehicle  2  coordinate  relative  to  burst  (kft) 

TSTAR 

i  nput 

RV  burst  time  relative  to  sea  level  impact 
(sec) 

Z1 

input 

Burst  elevation  (kft) 

IV 

input 

Vehicle  type 

S 

output 

Maximum  ratio  of  environment  level  to  blast 
kill  criteria 

IM 

output 

Number  of  maximum  environment 

C0MM0N  /HTRAJ/ 

input 

TIM,  V BAR,  ANGL,  MM 

C0MM0N  /WFRT/ 

output 

0PPK,  0DPK 

C0W0N  /FRAT/ 

input 

VI 

Routines  called: 

TERPL,  MATM62 ,  SCALKT ,  WFPK0P,  WFPK0D,  WFPKV, 
ORG 
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SUBROUTINE  LOAD  FLOW  DIAGRR  DIAGRAM 
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SUBR0UTINE  LP  ( XMW2 ,  GW) 


This  routine  calculates  the  local  properties  at  the  boundary  layer 
edge. 


XMW2 

input 

Mach  number  times  the  sine  squared  of  the 
angle  between  the  shock  and  RV  axis. 
(Presently  the  angle  is  90°) 

GW 

input 

Ratio  of  specific  heats 

behi nd 

the  shock 

C0MM0N 

/LPNC/ 

output 

All  variables 

C0MM0N 

/FSPC/ 

i  nput 

HT,  XM8,  P8PT8 

C0MM0N 

/  GPC/ 

input 

G,  GAM8,  XJ,  RGAS 

C0MM0N 

/SPC/ 

input 

PE,  PEPS 

I 


SUBROUTINE  MATM62  (TTY,  WSP,  CS,  WSR,  WST) 


This  routine  is  a  fit  to  the  1962  US  Standard  Atmosphere  and  was 
generated  by  the  Air  Force  Weapons  Laboratory. 


TTY 

input 

Altitude  (cm) 

WSP 

output 

Pressure  (dyne/ cm2) 

CS 

output 

Sound  speed  in  Air  (cm/sec) 

WSR 

output 

Density  (g/cm3) 

WST 

output 

Temperature  (°K) 

NZ 

data 

Number  of  data  values  in  the  atmosphere 
data  tables 

RH0Z 

data 

Ambient  sea  level  density  (g/cm3) 

TABAT(l) 

data 

Gas  constant  (8.3144xl07  ergs/ role/ deg) 

TABAT(2) 

data 

Radius  of  the  earth  (6.367488x10®  cm) 

TABAT(3) 

data 

Acceleration  of  gravity  at  sea  level 
(9. 80665x1 02  cm/ sec2) 

TABAT(4) 

data 

Molecular  weight  of  air  at  sea  level 
(28.9644) 

TABZC22) 

data 

Atmosphere  altitude  data 

TABL(21 ) 

data 

Atmosphere  molecular  scale  temperature 
gradient  (deg/cm) 

TABT(22) 

data 

Atmosphere  molecular  scale  temperature 
(°K/cm) 

TABP(22) 

da  ta 

Atmosphere  pressure  (dyne/cm2) 

C0NS 

data 

Combination  of  TABAT's 

C0NS  =  TABAT(3)*TABAT(4)*TA8AT(2)**2/TABAT(1 ) 

Routines  called: 

IBY 

SUBROUTINE  MESA G£  (MINOEX,  Ml,  N0,  RN0) 


Writes  error  and  key  diagnostic  messages 


MiNDEX 

i  nput 

F0RMAT  statement  number  for  the 
message 

desired 

Ml 

input 

Index  of  the  calling  subroutine 
in  array  LABEL  in  BL0CKDATA 

name  stored 

N0 

input 

Integer  variable  transferred  to 
message 

the  output 

RN0 

input 

Real  variable  transferred  to  the  output 
message 

C0MM0N  /DTH/ 

input 

LABEL 

SUBR0UTINE  MNCNCX  (ICNCX) 


This  is  a  sub-executive  routine  used  to  calculate  modified  Newtonian 
normal  and  axial  force  coefficients  for  composite  sphere-cone  vehicles 
as  a  function  of  angle  of  attack  with  vehicle  design  variables  used 
as  parameters. 


ICNGX 

output 

Flag  indicating  that  drag  coefficients  were 
input 

C0MM0N 

/BC/ 

input 

TCDEG,  SINTC,  C0STC,  TANTC, 

RN,  RB 

C0MM0N 

/PTMC/ 

input 

PI,  XRAD 

C0MM0N 

/DRAG/ 

input 

NM 

C0MM0N 

/DRAG/ 

output 

CN0,  CX0,  ALPHA,  XMACH,  NA, 

NM 

C0MM0N 

/ FSPC/ 

output 

XM8 

PMACH 

data 

Mach  numbers 

PALPHA 

da  ta 

Angle  of  attack  values 

IPM 

data 

Number  of  Mach  numbers 

I  PA 

data 

Number  of  angles  of  attack 

Routines  called: 

CDPCMP 
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SUBROUTINE  MRVSUM  (IPHASE,  IPASS,  IMRV,  RADMRV,  TMRV,  AMKV ) 


Creates  an  exclusion  region  encompassing  the  individual  burst 
contours  for  an  MRV  array.  Static  contours  have  a  circular  region 
for  boost  and  elliptical  for  reentry.  Blast  contours  are  merged 
into  the  rectangular  dynamic  region.  AMM  units  are  used  for  all 
vehicles. 


IPHASE 

input 

See  /RAD/ 

IPASS 

input 

Program  replication  index  for  MRV  option  3 

1  -  First  pass  —  first  burst  (downrange) 

2  -  Second  pass  —  second  burst  (uprange) 

IMRV 

input 

Type  of  MRV  pattern  modeled.  Three  different 
exclusion  region  models  are  defined  in  order 

to  minimize  region  sizes  for  different  pat¬ 
terns 

1  -  Three  or  more  bursts  having  identical 

yields  and  burst  times.  A  circular 
contour  is  generated 

2  -  Two  simultaneous,  equal -yield,  in-line, 

trajectory-plane  bursts.  A  rectangular 
contour  is  generated 

3  -  Two  different  yield,  different  altitude, 

in-line  bursts  with  the  uprange  bursf 
delayed  in  time.  Although  coding  fo- 
this  option  has  been  completed,  it  nu¬ 
merous  logic  paths  have  not  been  veri¬ 
fied  so  it  is  not  included  at  this  time. 
A  dummy  routine,  MRV3 ,  is  called  if 
MRV=3. 


RADMRV 

input 

Pattern  radius  for  IMRV=1  .  Burst  location  on 
the  patterns  does  not  matter.  One  half  the 
burst  separation  distance  for  IMRV=2  or  3. 

The  two  burst  cases  are  centered  about  the 
nominal  aimpoint  used  by  AMM  (NMi) 

TMRV 

input 

Time  separation  of  bursts  for  IMRV=3  (sec) 

AMRV 

input 

Altitude  of  the  second  burst  for  IMRV=3  (ft) 

C0MM0N  / BST/ 

input 

All 

C0MM0N  /BST/ 

output 

All  except,  NBL ,  NF0 

C0MM0N  / RVC/ 

input 

All 

C0MM0N  /RVC/ 

output 

All  except  NTM 

SUBROUTINE  MRVSUM  FLOW  DIAGRAM 


NRV  PATTERN  TYPE? 


IN-LINE  SURSTS  wrTH  DIFFERENT 
DS  AND  HOB'S,  SPACED  RAOHRV 

central  aimpcint 


BOOST  PHASE?  > 

' - 1 — - - ' 

NO 

_ 1 _ 

IF  HAVE  STATIC  DATA. 
EXPANO  REGION  SICE 


IF  HAVE  OYNAMIC  DATA, 
EXPANO  REGION  SIZE 


p  —  1 

CALL  MRV3 

DUMMY  ROUTINE  IN 
BASIC  NWEM  CODE 

BATA  PACKAGE 

I 

i 

t 

; _ 

CIRCULAR  PATTERN  WITH 
RADIUS,  RADMRV  _ 


REENTRY  PHASE? 


TWO  IDENTICAL.  IN-LINE 
BURSTS  SPACED  RADMRV 
FROM  CENTRAL  AIMPOINT 


YES 

I 

i 

IF  HAVE  STATIC  DATA, 
EXPAND  REGION  SIZE 


IF  HAVE  DYNAMIC  DATA. 
EXPAND  REGION  SIZE 


'  HAVE  BLAST  OR  ' 
^DUST/  EJECTA  DATA?  , 

YES 

_ » _ 

CALL  IN-LINE 

MERGE  BLAST  DATA 
INTO  DUST.'  EJECTA 
DYNAMIC  CONTOUR 


EXPANO  REGION  SIZE 


IF  HAVE  STATIC  ; 
EXPAND  REGION  RA, 


IF  HAVE  BLAS',  •' 

annular  res::1 


IF  HAVE  DUST/E JE 
EXPAND  DYNAMIC  RE 


SUBR0UTIME  NUCLD  (T,  W,  SH0B ,  HT,  HM,  RM.  RB,  RH0) 

This  routine  controls  the  dust  cloud  size  and  density  determination. 


T* 

i 

input 

Time  after  burst  (seel 

w 

input 

Burst  yield  (Mt ) 

Sh0B 

i  nput 

Burst  scaled  height-of-burst  (ft/ kt  3) 

HT 

output 

Dust  cloud  top  height  (kft) 

HM 

output 

Dust  cloud  middle  height  (kft! 

RM 

output 

Dust  cloud  middle  radius  (kft! 

RB 

output 

Dust  cloud  base  radius  (kft) 

RH0 

output 

Dust  cloud  density  (g.'cm3! 

C0MM0N  /C0NVR, 

input 

CVD 

C0MM0N  /THREAT/ 

i  nput 

N.  W0LD,  SH0LD ,  KDUST 

C0MM0N  /THREAT/ 

output 

I.  N 

C0MM0N  /STUFF/ 

input 

TER 

C0MM0N  / V0LLY/ 

Output 

V 

C0MM0N  /FARAM/ 

output 

WD 

Routines  cal  1 ed: 

WSTUFF,  EARLY.  LATE.  VOLUME ,  A  IRMAS 
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SUBROUTINE  NUCLD  FLOW  DIAGRAM 


SET  STORAGE 
INDEX 


CALL  WSTUFF 


COMPUTE  CASE 
CONSTANTS 


YES 


HOB  TOO 
GREAT  FOR 
.DUST  CLOUD?, 


YES 


NO 
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j  'SUBR0UTINE  0UTIN  (IPASS) 

i 

j  This  routine  outputs  various  levels  of  input  ana  associated  variable 

{  values  to  define  the  input  of  the  problem. 


IPASS 

input 

Flag  indicating  pass  number  for  MRV 
type 

third 

C0MM0N  /HTRNS/ 

input 

All  variables 

C0MM0N  / ITRNS/ 

input 

ITYPE,  NTRJ ,  IMRV,  INEUT,  IXRAY 

C0MM0N  /TRNS/ 

input 

BLST1 ,  BLST2,  BIST 3,  DUST!  ,  DUST2 , 

DUST3 , 

PEBB1,  PEBB2 ,  THER1  ,  THERF,  XNEUT1  ,  XGM1  , 
EMP1 ,  RVHT,  RVANG,  RVCL,  RVCYL,  RVFL,  YLD, 
XH0B2,  XNEUTX,  AFRAC,  GAMDS,  GAMPL,  RADMAX, 
RADTIM,  XTM1 ,  XM2 


1  C0MM0N  /DLAG/ 

input 

XCN ,  XCA,  XA1  ,  XTM2 ,  NA, 

NM 

i 

)  C0MM0N  /0UT/ 

input 

Q2 ,  Q3 ,  Q5 ,  Q6  ,  Q8,  Q9, 

IQ1,  IQ4,  IQ7 

|  C0MM0N  /TRAJ/ 

1 

input 

All  variables 

C0MM0N  /WPN/ 

input 

FN,  FX,  FTDG,  THERM 

C0MM0N  / XT ITL/ 

input 

LEV0UT 

I 

1 


SU8R0UTINE  0UTRIT  (JFLAG,  INO,  NDUM2) 


Edits  interim  and  final  output  from  the  radiation  and  exclusion 
region  contour  generation  subroutines  including  corunon  blocks 
/0UT/  ,  /RR1/  ,  / RR2/  ,  /STA/,  /  RRR/  ,  and  /WPN/. 


I  FLAG 

input 

Flag  denoting  the  block  of  output  desired 

1  -  C0MM0N  /RRR/ 

2  -  C0MM0N  /0UT/ 

3  -  C0MM0N  /RR1/ 

4  -  C0MM0N  / RR2/ 

5  -  C0MM0N  /SPA/ 

6  -  C0MM0N  /WPN/ 

7  -  DUMMY 

INO 

input 

Integer  variable  transfer  to  facilitate  edit 
or  include  in  edit 

NDUM2 

input 

A  second  integer  variable  to  facilitate  edit 
or  include  in  edit 

C0MM0N  /0UT/ 

input 

All  variables 

C0MM0N  /RR1  / 

i  nput 

All  variables 

C0MM0N  /RR2/ 

input 

All  variables 

C0MM0N  /RRR/ 

input 

All  variables 

C0MM0N  /STA/ 

input 

All  variables 

C0MM0N  /RAD/ 

i  nput 

All  variables 

C0MM0N  /DTH/ 

input 

NAME 

C0MM0N  /WPN/ 

input 

NAME 

KK 

data 

Flag  permiting  a  heading  for  all  detailed 
program  output  edit 

Routines  called: 

ERR0UT 
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SU8R0UTINE  PTINT  (H2,  V2,  G2,  R2,  T2,  CD,  Q8,  RVAM,  C0SG2,  AT,  VI,  G1  , 
R1 ,  T1 ,  D1 ,  El  ,  B1  ) 

This  routine  performs  the  individual  time  step  integration  of  the 
RV  trajectory. 


H2 

input 

Previous  RV  altitude  (ft) 

V2 

input 

Previous  RV  velocity  (ft/ sec) 

G2 

input 

Previous  RV  angle  to  the  local  horizontal 
(radians) 

R2 

input 

Previous  RV  ground  range  (ft) 

T2 

input 

Previous  time  (sec) 

CD 

input 

Drag  coefficient 

Q8 

input 

Dynamic  pressure  on  RV  during  time  step 
(slug/ft-sec2) 

RVAM 

input 

RV  reference  area  divided  by  the  mass 
( ft2/ si  ug) 

C0SG2 

output 

Cosine  of  the  RV  angle  to  the  local  horizontal 
C0S(G2) 

HI 

output 

Updated  RV  altitude  (ft) 

VI 

output 

Updated  RV  velocity  (ft/sec) 

G1 

output 

Updated  RV  angle  to  the  local  horizontal 
(radians) 

R1 

output 

Updated  RV  ground  range  (ft) 

T1 

output 

Update  time  (sec) 

D1 

output 

Acceleration  of  RV  along  flight  path  (ft/ sec2) 

El 

output 

Time  rate  of  change  of  RV  altitude  (ft/ sec) 

B1 

output 

Time  rate  of  change  of  RV  angle  to  the  local 
horizontal  (radian/sec) 

C0MM0N  /PTMC/ 

input 

XM0 ,  R0 

DELT 

data 

Maximum  integration  time  step  (sec) 

DVEL 

data 

Maximum  velocity  change  over  time  step 
interval  (ft/sec) 

SUBR0UTINE  PTLAST  (JJ,  KK,  RSTEP) 


Calculates  lethal  volume  size  along  a  path  directed  toward  a 
surface  point  corresponding  to  the  final  remaining  altitude  in 
either  the  upper  or  lower  half  of  the  volume.  SUBROUTINES  PT4 , 
PT17,  and  PTMAX  previously  calculated  volume  dimensions  at  the 
other  altitudes  (times). 


JJ 

input 

Environment  index 

1  -  Neutrons 

2  -  Peak  prompt  gamma  dose  rate 

3  -  X-rays 

4  -  Thermal  radiation 

KK 

output 

Provision  for  adding  a  flag  calling  the  error 
message  routine,  MESAGE,  from  RADIAT.  Not 
currently  used 

RSTEP 

input 

Range  step  size  used  to  find  the  range  to 
the  criterion  level  by  iterative  conver¬ 
gence  (km) 

C0MM0N  / RR1  / 

output 

T,  RH,  R,  RDEL,  A,  PERCNT 

C0MM0N  /RAD/ 

input 

CRIT,  ABURST,  AUPDN,  ALAST,  IMAX 

AP0INT 

Altitude  of  lethal  volume  surface  point  being 
defined.  Changes  as  location  estimate  is 
refined  during  the  iterations  (km) 

RNGE 

Range  to  the  lethal  volume  surface  point 
being  defined.  Changes  as  location  esti¬ 
mate  is  refined  during  the  iterations  (km) 

Routines  cal  led: 

RH0X,  FENV ,  TRAJEC 

See  PTMAX  flow  diagrams  for  similar  function  flow  of  PTLAST 
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SUBR0UT INE  PTMASS  (QH1 ,  QH1  ,  QV1 ,  QM,  QA,  QRN,  QRB,  QTC,  QRUFF,  QACYL , 
QAFLAR,  QHC0N,  QHCYL,  QHFLAR) 

This  routine  is  a  sub-executive  routine  which  directs  the 
calculation  of  the  RV  trajectory  data  and  drag  coefficients. 


QHl 

input 

RV  initial  reentry  altitude  (ft) 

QG1 

input 

RV  initial  reentry  angle  (deg) 

QV1 

input 

RV  initial  reentry  velocity  (ft/ sec) 

QM 

input 

RV  mass  (slugs) 

QA 

input 

RV  reference  area  (ft2) 

QRN 

input 

RV  nose  radius  (in) 

QRB 

input 

RV  base  radius  (in) 

QTC 

input 

RV  cone  half  angle  (deg) 

QRUFF 

input 

RV  surface  roughness  (in) 

QACYL 

input 

RV  cylinder  half  angle  (deg) 

QAFLAR 

input 

RV  flare  half  angle  (deg) 

QHC0N 

input 

RV  cone  length  along  axis  (in) 

QHCYL 

input 

RV  cylinder  length  along  axis  (in) 

QHFLAR 

input 

RV  flare  length  along  axis  (in) 

C0MM0N  /TRAJ/ 

input 

NSAVE 

C0MM0N  /TRAJ/ 

output 

TIME,  RANGE,  ALTX,  VEL,  FPA 

C0MM0N  /DR I DR/ 

output 

XTR,  ISTPAS 

C0MM0N  /DRAG/ 

input 

CX0,  XMACH,  NM 

C0MM0N  /DRAG/ 

output 

WN 

C0MM0N  /BC/ 

output 

All  variables 

C0MM0N  / FSPC/ 

i  nput 

XM8,  Q8 

C0MM0N  / FSPC/ 

output 

H2 ,  V2 

C0MM0N  /GPC/ 

output 

CPS,  G,  GAM8,  XJ,  PR,  RGAS 
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C0MM0N  /PTMC/ 
Routines  called: 


output  PI,  XRAD,  AM0,  R0 

MNCNCX,  ALTSET,  FSPX,  FSP,  DRAGX,  TERPL, 
PTINT,  TRAJET,  TRAST 


QACYL,  QAFLAR,  QHC0N,  QNCYL  AND  QHFLAR  are  not  presently  used  in  this  ver¬ 
sion  of  the  NWEM  code 


SUBROUTINE  PTMASS  FLOW  DIAGRAM 


TES 


_ CALL  TRAST _ 

INVERT  TI*  AMO  RANGE  VECTORS 
TO  K  DEFINED  IN  TERMS  OF 
TltC  AND  RANGE  TO  IMPACT 

C  ReIirn  ) 
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5UBR0UT I NE  PTMAX  (JJ,  KK,  ASTART,  RSTEP,  UMAX) 


Calculates  lethal  volume  size  for  the  altitude  at  which  the 
trajectory  grazes  the  surface,  either  above  or  below  the  burst 
altitude.  Follows  calculation  of  top,  bottom,  and  burst  altitude 
dimensions  by  SUBR0UTINES  PT4  and  PT17.  Called  by  RADIAT  except 
for  a  surface  burst  or  if  the  grazing  trajectory  below  the  burst 
point  were  located  below  the  earth's  surface. 


input 


ASTART 


RSTEP 


ITMAX 


C0MM0N  / RR1  / 
C0MM0N  /RAD/ 
C0MM0N  /RAD/ 


Environment  index 

1  -  Neutrons 

2  -  Peak  prompt  gamma  dose  rate 

3  -  X-rays 

4  -  Thermal  radiation 

Provision  for  adding  a  flag  calling  the  error 
message  routine,  MESAGE,  from  RADIAT.  Not 
currently  used 

Initial  estimate  of  grazing  trajectory  con¬ 
tact  point  altitude  based  on  burst  and 
either  top  or  bottom  altitudes  of  the  lethal 
vol  ume 

Range  step  size  used  to  find  the  criterion 
level  by  iterative  convergence  (km) 

Maximum  number  of  time  steps,  i.e.,  the 
index  of  the  bottom  of  the  lethal  volume. 

Used  to  define  the  time  index  of  the  point 
determined  for  a  grazing  trajectory  below 
the  burst  point 


output  T,  RH,  R,  RDEL,  A,  PERCNT 
input  CRIT,  ABURST,  PI,  AUPDN,  IMAX 
output  ALAST 


Altitude  of  lethal  volume  surface  point 
being  defined.  Changes  as  location  estimate 
is  refined  during  the  iteration  (km) 

Range  to  the  lethal  volume  surface  point 
being  defined.  Changes  as  location  estimate 
is  refined  during  the  iterations  (km) 


output 


input 


input 


input 


I MpWMft*' 


ELEV 


Routines  cal  led: 


input  Elevation  angle  of  the  trajectory  w.r.t.  the 
local  horizontal.  Changes  as  location  esti¬ 
mate  is  refined  during  the  iterations  (radians; 

PH0X,  FENV,  TRAJEC 


SUBROUTINE  PTMAX  FLOW  DIAGRAM  * 


*  Functional  flow  for  FTU,  FT«  ond  FTtAST  H  »1«11tr. 
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SUBROUTINE  PT17  (JJ,  KK,  ITMAX) 


Calculates  lethal  volume  size  along  the  ground  zero  (vertical  axis) 
in  either  direction.  If  the  surface  is  encountered  before  the 
downward  size  is  determined,  the  range  to  a  point  on  the  earth's 
surface  is  determined  neglecting  radiation  scattering  effects  from 
the  surface  itself. 


JJ 


KK 


ITMAX 


input  Environment  index 

1  -  Neutrons 

2  -  Peak  prompt  gamma  dose  rate 

3  -  X-rays 

4  -  Thermal  radiation 

output  Flag  directing  that  the  calculation  be 

terminated  after  MESAGE  is  called  to  print  a 
message  that  the  problem  does  not  converge 

output  Index  of  final  time  point  (see  C0MPAR) 


C0MM0N  / RR1 /  output  T,  RH,  R,  RDEL ,  A,  PERCNT 

C0MM0N  /RAD/  input  CRIT,  ABURST,  AUPDN,  ICHK 


AP0INT 

RNGE 

STEPA 


RSTEP 


Altitude  of  lethal  volume  surface  point  being 
defined.  Changes  as  location  estimate  is 
refined  during  the  iterations  (km) 

Range  to  the  lethal  volume  surface  point  being 
defined.  Changes  as  location  estimate  is  re¬ 
fined  during  the  iterations  (km) 

Altitude-dependent  estimate  of  range  step 
size  sufficiently  small  that  an  iterative 
convergence  to  the  correct  lethal  volume 
size  is  possible  (km) 

Range  step  size  used  to  find  the  exact  range 
to  the  criterion  level  by  iterative  conver¬ 
gence  (km) 


Routines  called: 


RH0X,  FENV,  TRAJEC,  MESAGE 


See  PTMAX  flow  diagrams  for  similar  function  flow  of  PT1 7 
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SUBR0UTINE  PT4  (JJ,  KK) 


Calculates  the  lethal  volume  size  at  the  burst  altitude  and  assiqns 
the  descriptive  data  to  index  4  in  the  time  array. 


JJ 


KK 


input  Environment  index 

1  -  Neutrons 

2  -  Peak  garrnia  dose  rate 

3  -  X-rays 

4  -  Thermal  radiation 

output  Flag  directing  that  RADIAT  call  MESAGE  to 
print  an  error  message  indicating  that  the 
problem  cannot  converge  to  the  criterion  en¬ 
vironment  level.  It  then  terminates  the 
radiation  environment  calculation 


C0MM0N  / RR1 /  output  T,  RH,  R,  RDEL,  A,  PERCNT 

C0MM0N  /RAD  input  CRIT,  ABURST,  ICHK 


RSTEP 


RNGE 


Routines  called: 


Range  step  size  used  to  find  the  exact  range 
to  the  criterion  level  by  iterative  conver¬ 
gence.  An  estimate  of  one-half  the  exoatmos 
pheric  range  to  the  criterion  level  is  used 
initially  for  burst  altitudes  above  20  km 
and  an  altitude-varying  size  at  lower  alti¬ 
tudes  (km) 

Range  to  the  lethal  volume  surface  point 
being  defined.  Changes  as  location  estimate 
is  refined  during  the  iterations  (km) 

RH0X,  FENV,  TRAJEC 


See  PTMAX  flow  diagrams  for  similar  function  flow  of  PT4 
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SUBROUTINE  QUAD  (X,  XT,  YT,  Y) 


This  is  a  quadratic  interpolation  routine  for  three  point 
arrays.  No  search  is  performed,  nor  is  extrapolation  done 
points  values  are  returned. 


X 

input 

XT 

input 

YT 

input 

Y 

output 

Independent  variable 
Independent  variable  array 
Dependent  variable  array 
Dependent  variable  value 
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pai r  wise 
.  End 


SUBROUTINE  RADIAT  (IP,  I PASS,  AB,  CC1  ,  CC2 ,  CC3,  CC4) 


Sab-executive  routine  called  by  the  main  program  to  process  radiation 
environments.  Calculates  lethal  volumes  and  exclusion  regions  at  up 
to  7  times  for  each  radiation  environment.  Currently  evaluates  four 
environments:  prompt  neutron  fluence,  prompt  gamma  peak  dose  rate, 
X-ray  fluence,  and  thermal  radiation  fluence.  A  single  encompassing 
exclusion  region  is  returned  for  either  boost  or  reentry  phase. 


IP 

input 

Flight  phase  index,  equal  to  IPHASE 

2  -  Boost  phase 

4  -  Reentry  phase 

IPASS 

input 

Flag  indicating  which  burst  is  being  calculated 
for  MRV  Option  3 

0  -  First  burst  (downrange  w.r.t.  RV  launch 
point) 

1  -  Second  burst  (uprange  w.r.t.  RV  launch 
point) 

AB 

input 

Burst  al titude  (ft) 

CC1 

input 

Environment  criterion  for  neutrons,  equal  to 
CRIT(l)  (neutrons/cm2) 

CC2 

input 

Environment  criterion  for  prompt  gammas,  equal 
to  CRIT(2)  ( rad ( Si )/sec) 

CC3 

input 

Environment  criterion  for  X-rays,  equal  to 
CRIT(3)  (cal /cm2) 

CC4 

input 

Environment  criterion  for  thermal  radiation, 
equal  to  CRIT (4)  (cal/cm2) 

C0MM0N  /RRR/ 

output 

UMAX,  TIMIN,  JINDEX,  MAXIN,  IL00PS,  ITM 

C0MM0N  / RR1  / 

Output 

R,  A 

C0MM0N  /RAD/ 

output 

CRIT,  ABURST,  PI,  AUPDN,  ISKIP,  IPHASE,  ICHK 

C0MM0N  /RAD / 

i  n  pu  t 

IMAX 

C0MM0N  / XTITL/ 

input 

LEV0UT 

STEPA 

Estimated  step  size  defined  such  that  the 
scheme  for  environment  definition  will  converge 

Routines  called: 

SOURCE,  0UTRIT,  MESAGE,  PT4,  PT17,  PTLAST, 

PTMAX,  C0MPAR,  RADMAX 
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SUBROUTINE  RADMAX  (IPHASE,  NENV,  ISKIP) 


Determines  the  maximum  surface  exclusion  region  resulting  from  all 
radiation  environments  over  all  times. 


IPHASE 


ISKIP 


C0MM0N  / DTH 
COMMON  /XT TL/ 
COMMON  / RRR/ 
COMMON  / RR1 / 
C0MM0N  / RR2/ 
COMMON  /STA/ 


input 


input 


input 


input 


input 


Routines  called: 


FI ight  phase  index 
2  -  Boost  phase 
4  -  Reentry 

Total  number  of  radiation  environments  — 
currently  4 

Flag  for  each  environment,  set  equal  to  1  if 
the  environment  is  to  be  shipped  (environment 
index  order  as  for  CR IT ) 


input  LEV0UT 
input  All  variables 

input  All  except  PERCNT 

output  All  variables 

output  All  variables 


Time  index  plus  fractional  increment  towards 
the  next  index  that  the  time  corresponding 
to  the  maximum  horizontal  radius  (over  all 
environments)  occupies  in  the  final  time 
array.  Used  to  adjust  one  of  the  final 
time  array  points-,  TEX(I),  to  coincide  with 
the  time  of  maximum  horizontal  radius 

Linear  fraction  of  a  time  interval  T(I)-T(I+1 
that  TEX(J)  occupies.  Defined  by  SUBROUTINE 
TERP  such  that  linear  interpolation  of  the 
exclusion  contour  data  for  the  environment 
in  question  can  be  performed  to  compare  it 
with  data  from  other  environments  at  the 
same  final  time  so  that  the  maximum  exclusion 
contour  parameters  can  be  obtained 


TERP,  RAD2 ,  0UTRIT 
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SUBROUTINE  RADMAX  FLOW  DIAGRAM 


OATA  TOR  A  SINGLE  ENVIRONMENT? 

i 

NO 


TRANSFER  CONTOUR  OATA 
TO  FINAL  ARRAY 


MAXIMUM  HORIZONTAL  CONTOUR  SIZE 
AT  BOTTOM  OF  VOLUME  7 
(LIKELY  FOR  SURFACT  BURST) 


ywximum  horizontal  contour  size 

AT  TOP  OF  VOLUME  ?  (UNLIKELY) 


OIVIOE  CONTOUR  SIZE  ARRAY 
INTO  EQUAL  TIME 
INCREMENTS 


SELECT  7  EQUISPACED  TIMES  SPAWNING 
EXCLUSION  CONTOUR  DURATION 


ADJUST  THE  TIME  STEP  CLOSEST  TO  The 
MAXIMUM  CONTOUR  HORIZONTAL  EXTENT 
TO  CONINCIUE  WITH  THAT  r»Mf 


LOOP  OVER  ALL  RADIATION  ENVIRONMENTS 


I  LOOP  OVER  THE  ADJUSTED  TIME  ARRAT 


VALIO  OATA  FOR  THIS  ENVIRONMENT’ 

■  I 

YES 


CURRENT  ADJUSTED  TIME  LIE  WITHIN 
TIME  ARRAT  FOR  THIS  ENVIRONMENT’ 


CONTOUR  PARAMETERS  PREVIOUSLY 
DEFINED  FOR  THIS 
ADJUSTED  TIME? 


t 

DEFINE  EXCLUSION  CONTOUR 
SIZE  ANO  TRANSLATION  FOR 

1 

CALCULATE  MAXIMUM  HORIZONTAL 
RADIUS  AT  THIS  ADJUSTEl.  TIME j 

REPUCE  STORED  CONTOUR  PARAMETER 
DATA  WITH  OATA  FOR 
_ CURRENT  ENVIRONMENT _ 


RADIUS  EXCEED  PREVIOUSLY 
STORED  RADIUS’ 


CALL  RAD2 


SELECT  EXCLUSION  CONTOUR 
MAXIMUM  DIMENSIONS 


LMD  COWON  BLOCK  FOR  STATIC 
CONTOUR  RADIATION /BLAST  MERGING 


SUBR0UTINE  RAD2  (NAXITM,  IPHASE) 


Selects  the  maximum  values  of  radiation  exclusion  region  data, 
identifying  each  causative  environment.  Converts  exclusion  re¬ 
gion  times  to  negative  times. 


MAXITM 


IPHASE 


C0MM0N  / RR2/ 
C0MM0N  / RR2/ 


input  Maximum  number  of  times  specified.  Redefined 
as  7  in  RADMAX  if  2  or  more  radiation  environ¬ 
ments  have  been  calculated 

input  Flight  phase 

2  -  Boost 

4  -  Reentry 

input  TEX,  RHEX,  RUPEX,  RDNEX,  NAMEX 

output  TEX,  TMIN,  TMAX,  RHMAX ,  RUPMIN,  RDNMAX , 

NAMET1 ,  NAMET2 ,  NAMEH ,  NAMEU,  NAMED 


FUNCTI0N  RBASE  (T) 

This  routine  defines  the  dust  cloud  base  radius  after  fireball 
breakaway  (kft). 


T 

input 

Time 

after  burst  (sec) 

C0MM0N 

/THREAT/ 

input 

I 

C0MM0N 

/STUFF/ 

input 

ABO, 

BBO,  RBC,  RBS,  TBC,  TBS 

H 


k 


t. 
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SUBROUTINE  REENTR  (T,  X2 ,  Y2,  22,  20,  Z1  ,  YDT ,  YSTAR,  TSTAR ,  S,  IV,  IM) 


This  routine  controls  search  for  RV  intercepts wi  th  shock  waves  and 
defines  the  resulting  environment  severity. 


T 

input 

X2 

i  nput 

Y2 

input 

12 

i  nput 

10 

input 

11 

input 

YDT 

i  nput 

YSTAR 

i  nput 

TSTAR 

input 

s 

output 

r  \  ; 

*  V 

input 

IM 

output 

C0MM0N  /HTRAJ/ 

input 

Routines  called: 

Time  of  RV2  impact  after  RV1  burst  (sec) 


Position  of  RV2  burst  relative  to  RV1  burst, 
x  component  (kft) 


Position  of  RV2  burst  relative  to  RV1  burst, 
y  component  (kft) 


RV2  burst  elevation  (kft) 
Ground  plane  elevation  (kft) 
RV1  burst  elevation  (kft) 
Effective  blast  yield  (fit) 


Y  offset  of  RV2  due  to  burst  elevation  above 
sea  level  (normal  trajectory  termination) 
(kft) 


Time  offset  (along  trajectory)  of  RV2  due  to 
burst  elevation  above  sea  level  (normal 
trajectory  termination)  (sec) 


Environment  severity.  Ratio  of  intercepted 
environment  to  kill  criteria. 


Vehicle  type:  2  for  booster,  4  for  RV 


Number  of  blast  environment  with  the  maximum 
severi ty . 


TBP,  RBP ,  ABP,  NT RAJ 
TERPL,  RSHK,  VL0C ,  L0AD 
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SUBROUTINE  REENTR  FLOW  DIAGRAM 


l 

t 
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SUBROUTINE  REXBDY  (XM,  REX2RN) 

This  routine  calculates  the  local  Reynold's  number  at  a  position 
along  the  RV  axis  two  times  the  RV  nose  radius. 


XM 

input 

Mach  number 

REX2RN 

output 

Local  Reynolds  number  at  a  position  along 
the  RV  axis  twice  the  RV  nose  radius 

C0MM0N  /BC / 

input 

TCRAD,  SINTC,  C0STC ,  TANTC,  RN 

C0MM0N  /FSPC/ 

input 

P8,  HT 

C0MM0N  / GPC/ 

input 

G,  GAMS,  XJ 

C0MM0N  /PTMC/ 

input 

PI 

GAM3 

data 

Coefficient  in  correlation  at  the  elevated 
temperature  (.718) 

FUNCTION  RH0X  (Z1 

,  Z2 ,  R) 

Determines  the  air  mass  integral  in  gm/cm2  between 
Based  upon  flat  earth  geometry.  If  the  two  points 
MATM62  is  called  to  return  air  density. 

two  points, 
are  coal titude, 

Z1 

input 

Altitude  of  first  point  (km) 

Z2 

input 

Altitude  of  second  point, 
in  any  order  (km) 

Z1  and  Z2  can 

be 

R 

input 

Slant  range  between  Z1  and 

Z2  (km) 

Z  (92 ) 

data 

Altitude  table  (kft) 

RH0DZ(92 ) 

data 

Vertical  air  mass  integral 
corresponding  altitude,  Z, 

from  0  kft  to 
(gm/cm2) 

the 

Routines  called: 

TERPL,  MATM62 
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FUNCTI0N  RSHK  (TIM,  22  ,  Zl,  VLD) 


This  routine  calculates  the  shock  radius  (kft)  from  the  burst  point. 
It  is  based  on  the  AFWL  1  kt  standard. 


TIM 

input 

Time  after  burst  (sec) 

12 

input 

Elevation  of  the  point  of  interest  (kft 

Zl 

input 

Burst  elevation  (DUMMY) 

YLD 

input 

Burst  yield  (Mt) 

B 

data 

Constant  for  curve  fit 

C 

data 

Constant  for  curve  fit 

CZ 

data 

Constant  for  curve  fit 

BZ 

data 

Constant  for  curve  fit 

TA3 

data 

Constant  for  curve  fit 

TPWR1 

data 

Constant  for  curve  fit 

TPWR2 

Routines  called: 

data 

Constant  for  curve  fit 

SCALKT 

SU8R0UTINE  SCALKT  (H,  HFPT,  WE,  VSCALE,  DSCALE ,  TSCALE,  CSCALE,  PSCALE) 

This  routine  sets  up  scale  factors  for  modified  Sachs  scaling  of 

1  kT  sea  level  data. 


H 

i  nput 

Burst  elevation  (dummy) 

HFPT 

i  nput 

Elevation  of  point  of  interest  (cm) 

WB 

input 

Burst  yield  (kt) 

VSCALE 

output 

Velocity  scale  factor 

DSCALE 

output 

Density  scale  factor 

TSCALE 

output 

Time  scale  factor 

CSCALE 

output 

Length  scale  factor 

PSCALE 

output 

Pressure  scale  factor 

C0MM0N  /C0NST/ 

input 

THRO,  PI,  Cl,  R1 

Routines  called: 

MATM62 

SUBR0UTINE  SIZER  (EM) 

This  routine  sizes  the  grid  for  subsequent  definition  of  intercepted 

kinetic  energy  points. 

EM  output 

C0MM0N  /BURST/  input 

C0MM0N  /DIMEN/  output 

C0MM0N  /HTRAJ/  input 


Trajectory  slope  at  cloud  top 

HT,  HM,  RM,  HBRV 
NY,  NX,  YMAX,  DY ,  DX 
AA,  RR,  NN 


Routines  called: 


SL0PE ,  TERPL 


FUNCTION  SL0PE  (Z) 


This  routine  defines  the  vehicle  trajectory  slope  at  a  given 

al  titude. 


Z 


input  Altitude  (kft) 


C0MM0N  /HTRAJ/ 


input  AA,  RR,  NN 


SUBROUTINE  S0L  V  ^  TC ,  XM,  BR,  C,  N,  ANS) 

this  routine  performs  the  third  and  fourth  order  interpolations 
required  by  the  correlation  function  of  the  drag  coefficient. 


TC 

input 

RV  cone  hal f  angle  (deg) 

XM 

i  nput 

Mach  number 

BR 

input 

RV  bluntness  ratio  or  log  of 
ratio 

the  bluntness 

C 

input 

Correlation  data 

N 

input 

Order  of  the  correlation  fit 
angle 

i n  RV  cone  hal f 

ANS 

output 

Correlation  function  result 
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SUBROUTINE  S0URCE 


Generates  weapon  radiation  output  source  strengths  using  data  from 
/WPN/  to  load  /0UT/  so  that  the  radiation  environment  at  a  point 
can  be  calculated  by  FENV.  Converts  air  transmission  data  in  TX, 
TN,  and  TG  arrays  to  1 og, Q  so  that  semi  log  interpolation  can  be 
used. 


C0MM0N  /0UT/ 

input 

TX,  TN,  TG,  N0X,  N0N ,  N0G 

C0MM0N  /0UT/ 

output 

SX,  TX,  SN,  TN,  SG,  TG,  N0X ,  N0N,  N0G 

C0MM0N  /RAD/ 

input 

ISKIP 

C0MM0N  /RAD/ 

output 

W0UT 

C0MM0N  /WPN/ 

i  nput 

All  variables 

C0MM0N  /WPN/ 

output 

FN,  FG,  FX,  FT,  DELTG 

SUBR0UTINE  SP 

This  routine  calculates  the  static  pressure  and  static  pressure  to 
ambient  pressure  ratio  for  the  RV  drag  coefficient  calculation. 


C0MM0N  /BC/ 

input 

SINTC 

C0MM0N  / FSPC/ 

input 

P8 ,  XM8 

C0MM0N  /GPC/ 

input 

GAM8 

C0MM0N  /SPC/ 

output 

All  variables 
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SUBROUTINE  STATIC 


Combines  the  overall  radiation  static  contour  with  the  negative 
time  (i.e.,  launch  before  burst)  portion  of  the  blast  contour  to 
get  a  net  static  region  exclusion  contour  for  the  boost  pnase. 
Dust  and  ejecta  are  assumed  non-limiting  at  negative  times.  If 
the  exclusion  region  does  not  encompass  ground  zero,  it  is  as¬ 
sumed  that  increasing  the  size  of  the  region  to  include  all  the 
projected  area  would  be  excessive  since  an  offset  region  is  pro¬ 
hibited  by  AMM,  so  a  message  is  printed  and  no  static  contour  is 
produced.  This  limits  the  current  modeling  to  low  altitude 
bursts  whose  lethal  volumes  reach  the  surface. 


C0MM0N  / BST/ 
C0MM0N  / DTH/ 
C0MM0N  /ST A/ 
C0MM0N  /'XT ITU- 


output  F0RAUI,  F0TIMS 

input  NAME 

input  All  variables 

input  LEV0UT 


I  PASS 


RMAX 


NRMAX 

RTEST 


Flag  indicating  environments  being  processed 
or  that  have  been  processed,  either  1  or  2 
at  this  time.  Note,  this  is  different  from 
IPASS  as  used  in  the  main  program 

output  Maximum  distance  (positive)  from  ground  zero 
to  any  point  on  the  exclusion  contours, 
equal  to  F0RADI  when  converted  to  nmi  if 
ground  zero  is  encompassed  (km) 

Label  of  limiting  environment  (radiation  or 
blast)  (bollerith  variable) 

output  If  ground  zero  is  not  encompassed,  the  closest 
approach  of  the  contour  to  ground  zero  is  out¬ 
put  by  MESAGE  (km) 


Routines  called: 


MESAGE,  0UTRIT 


SUBROUTINE  STATIC  FLOW  DIAGRAM 
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SUBR0UTINE  STATRV 


Converts  the  static  exclusion  contour  parameters  for  the  radiation 
environment  (only  static  region  contributor)  to  AMM  output  format 
The  duration  of  the  region  is  overestimated  if  the  lethal  volume 
doesn’t  intersect  the  surface  since  time  offset  is  not  permitted 


C0MM0N  / RVC/ 
C0MM0N  /STA/ 
C0MM0N  mni/ 


output  FRSLEN,  FRSWID,  FRSCEN,  TIMST 
input  All 
input  LEV0UT 


Routines  called: 


0UTRIT 


FUNCTI0N  TBREAK  (W) 

This  routine  provides  a  characteristic  time  (min)  used  in  definition 
of  the  dust  cloud  base  radius. 

input  Burst  yield  (Mt) 

data  Minimum  yield  ( .0353) 

data  Maximum  yiel d  (200.0) 

^  data  Yield  array 

B  data  Characteristic  time  array 
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SUBROUTINE  TERP  (TGIVN,  III,  J,  I  FIND ,  FRAX) 


Finds  the  index  of  the  lower  bound  in  the  T  array  interval  that 
TGIVN  is  within  plus  the  linear  fraction  of  the  interval  to 
TGIVN.  T  is  monotonical 1y  decreasing  and  TGIVN  lies  within  the 
T  array.  The  index,  IFIND,  and  fraction,  FRAX,  are  used  to 
interpolate  region  parameters  which  vary  with  T. 


TGIVN 

input 

Time  between  lethal  volume  point  and  the 
surface  which  lies  within  the  time  array, 

T ( I , J )  (sec) 

III 

input 

Maximum  time  index,  I,  of  the  time  array, 

T( I , J ) 

J 

input 

Radiation  environment  index 

1  -  Neutrons 

2  -  Prompt  gamma 

3  -  X-rays 

4  -  Thermal  radiation 

IFIND 

output 

Lower  bound  index  of  the  time  array  interval 
containing  TGIVN 

FRAX 

output 

Linear  fraction  of  the  interval 

T(I+1,0)  -  T(I,J)  corresponding  to  the  time 
between  T(I,J)  and  TGIVN 

C0MM0N  /  RR1  / 

input 

T 
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FUNCTI0N  TERPL  (UX,  X,  Y,  N,  NF1,  NF2) 


This  function  is  a  linear  interpolation  routine  for  pair  wise  arrays 
For  a  given  set  of  array  data  and  a  value  somewhere  in  the  range  of 
the  independent  variable,  it  returns  the  associated  dependent  value 
by  linearly  interpolating  in  the  data  arrays.  For  a  value  outside 
the  independent  variable's  range  the  function  returns  a  linearly 
extrapolated  value  or  the  end  point  value. 


UX 

input 

Input  value  of  the  independent  variable 

X 

input 

Independent  variable  array  in  ascending 
order 

Y 

input 

Dependent  variable  array 

N 

input 

Number  of  data  values 

N  FI 

input 

Interpolation  flag  for  index  value  in 
independent  array  (C/1  -  unknown/known) 

NF2 

i  nput 

Extrapolation  flag  (0/1  -  linear  extrapola 
ti on/ end  value) 

Routine  called: 

IBY 

FUNCTION  TERPL2  (UX,  UY,  X,  Y,  Z,  NX,  NY,  NF) 


This  function  is  a  linear  interpolation  routine  for  two  dimensional 
interpolation  in  an  array.  For  a  given  set  of  array  data  and  set 
of  values  somewhere  in  the  ranges  of  the  independent  variables,  it 
returns  the  associated  dependent  value  by  linearly  interpolating 
in  the  data  arrays.  For  values  outside  the  ranges  of  the  indepen¬ 
dent  variables  the  function  returns  a  linearly  extrapolated  value 
or  the  end  point  value. 


UX 

input 

Input  value  of  first  independent  variable 

UY 

input 

Input  value  of  second  independent  variable 

X 

input 

First  independent  variably  array  in  ascen¬ 
ding  order 

Y 

input 

Second  independent  variable  array  in  ascen¬ 
ding  order 

Z 

input 

Two  dimension  dependent  variably  array 

NX 

input 

Number  of  data  values  of  first  independent 
variable 

NY 

input 

Number  of  data  values  of  second  independent 
variable 

NF 

input 

Extrapolation  flag  (0/1  -  linear  extrapolation/ 
end  values) 

I  BY 


Routines  called: 


SUBR0UTINE  TNFR  (N,  X,  V) 


This  routine  transfers  data  from  one  array  to  another  array. 

input  Number  of  data  to  be  transferred 
input  Data  array  to  be  transferred  from 
output  Data  array  to  be  transferred  to 


SUBR0UTINE  T0NSET 

This  routine  determines  the  altitude  at  which  turbulence  comes  onto 
the  RV  body  for  the  drag  calculation. 


C0MM0N  /DR I DR/ 

output 

IT  J  F 

C0MM0N  /BC/ 

input 

TCDEG,  TCRAD,  SINTC,  C0STC,  RN,  R6 

C0MM0N  / FSPC/ 

This  conroon  is  used  in  a  temporary  data  pas¬ 
sing  mode  during  iterations  to  determine  tur 
bulent  onset  attitude 

C0MM0N  / LPNC/ 

input 

XME2,  REFT2 

C0MM0N  /GPC / 

input 

GAM8,  RGAS 

C0MM0N  /TRNSC/ 

input 

LAMINR 

C0MM0N  /TRNSC/ 

output 

LAM IN R,  XKTR 

C0MM0N  ,/PTMC/ 

input 

PI 

GAM1 

data 

Coefficient  in  correlation  equation 
local  Reynolds  number  (901000.) 

for 

GAM2 

data 

Coefficient  in  correlation  equation 
local  Reynolds  number  (-0.784) 

for 

12 

data 

Constant  2  used  in  iteration  scheme 

Routines  called: 

ERR0UT,  REXBDY,  TPPTID,  TERPL,  RSP, 
FNFRE 

FNFME, 
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SUBROUTINE  TRAJEC  (IPHASE,  M0OE,  AAA,  RD,  TT,  ELEV,  VEL) 


Returns  trajectory  data  for  either  boost  or  reentry  phase  after 
entering  with  altitude  or  time  interval  to  the  surface.  Data 
generated  consist  of  time  to  surface  (positive),  ground  range 
along  the  earth's  surface,  altitude,  velocity,  and  angle  with 
respect  to  the  horizontal.  A  warning  message  is  written  if  the 
time  or  altitude  is  outside  the  array  bounds  and  extrapolation 
is  used  to  obtain  the  data. 


IPHASE 


input 


Flight  phase  -  not  currently  used  since 
trajectory  data  are  stored  and  used  in  an 
identical  manner  for  each  phase 
2  -  Boost 
4  -  Reentry 


M0DE 


input 


Flag  indicating  entry  variable 

1  -  Enter  with  altitude 

2  -  Enter  with  time  to  surface 


AAA 

RD 


input  or 
output 


output 


TT 

ELEV 


input  or 
output 


output 


VEL 

ALT 


output 

output 


Routines  called: 


Altitude  (km) 


Ground  range  between  trajectory  point  and 
the  surface.  Approximately  equal  to  a 
horizontal  range  for  the  distances  involved, 
(km) 


Flight  time  between  surface  and  the  trajec¬ 
tory  point  (sec) 


Elevation  angle  with  respect  to  the  local 
horizontal  (radians) 


Vehicle  velocity  (km/sec) 


Altitude,  AAA,  (with  units  changed)  printed 
by  MESAGE  if  the  altitude  is  outside  the 
/TRAJ/  array  bounds  (ft) 


MESAGE,  TERPL 
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SU8R0UTINE  TRASET  (H2,  ALT,  D1,  El,  B1 ,  T2,  V2,  G2,  C0SG2 ,  R2,  TIME, 
VEL,  ANGL,  RANGE) 


This  routine  performs  the  time  step  integration  of  the  RV  trajectory 
at  the  altitude  break  points. 


H2 

input 

Previous  RV  altitude  (ft) 

ALT 

input 

Altitude  break  point  (ft) 

D1 

input 

Acceleration  of  RV  along  flight  path 
(ft/ sec2 ) 

El 

input 

Time  rate  of  change  of  RV  altitude  (ft/sec) 

B1 

input 

Time  rate  of  change  of  RV  angle  to  the  local 
horizontal  (radian/ sec) 

T2 

input 

Previous  time  (sec) 

V2 

input 

Previous  RV  velocity  (ft/sec) 

G2 

input 

Previous  RV  angle  to  the  local  horizontal 
(radian) 

C0SG2 

input 

Cosine  of  the  RV  angle  to  the  local  horizontal 
C0S(G2) 

R2 

input 

Previous  RV  ground  range  (ft) 

TIME 

output 

Time  at  altitude  break  point  (sec) 

VEL 

output 

RV  velocity  at  altitude  break  point  (ft/ sec) 

ANGL 

output 

RV  angle  to  the  local  horizontal  at  altitude 
break  point  (deg) 

RANGE 

output 

RV  ground  range  at  altitude  break  point  (NMi ) 

C0MM0N  / PTMC/ 

input 

XRAD,  R0 
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SUBR0UTINE  TRAST  (TIME,  RANGE,  NSAVE) 

This  routine  inverts  the  order  and  value  of  the  calculated  time  and 
range  arrays  of  the  R V  trajectory. 


TIME 

input 

RV  trajectory  time  data  (sec) 

output 

RANGE 

input 

RV  trajectory  ground  range  (NMi) 

output 

NSAVE 

input 

Number  of  RV  trajectory  data 

SUBR0UTINE  TRPTID 


This  routine  computes 
the  drag  calculation 

the  laminate  turbulent  transition  point  for 
in  the  non-coupled  mode  and  through  T0NSET. 

C0MM0N  / DRI DR/ 

output 

XTR 

C0MM0N  /BC/ 

input 

TCDEG,  TCRAD,  SINTC,  C0STC,  RN  RB 

C0MM0N  /FSPC/ 

input 

XM8 

C0MM0N  /LPNC/ 

input 

XME2,  REFT2 

C0MM0N  /TRNSC/ 

output 

LAM 1 NR,  RNRTR 

C0MM0N  /PTMC/ 

input 

PI 

REXMIN 

data 

Minimum  Reynolds  number  allowed  for  tran¬ 
sition  (5x10s) 

XKTR 

data 

Transition  point  parameter  constant  (41043.) 

Routines  called: 


FNFME,  FNFRE 


SUBROUTINE  VL0C  (TB 

,  X,  Y, 

2,  X2,  Y2,  Z2  ,  YSTAR,  TSTAR,  IV) 

This  routine 

locates 

the  vehicle  at  the  given  time. 

TB 

input 

Time  to  burst  pt  (RV)  or  time  since  launch 
(booster)  (sec) 

X 

output 

Vehicle  location,  x  coordinate  (kft) 

Y 

output 

Vehicle  location,  y  coordinate  (kft) 

2 

output 

Vehicle  location,  z  coordinate  (kft) 

X2 

input 

Burst  or  launch  point,  x  coordinate  (kft) 

Y2 

input 

Burst  or  launch  point,  y  coordinate  (kft) 

72 

input 

Burst  or  launch  point,  z  coordinate  (kft) 

YSTAR 

input 

Trajectory  range  offset  due  to  RV  burst  above 
sea  level  (kft) 

TSTAR 

input 

Trajectory  time  offset  due  to  RV  burst  above 
sea  level  (kft) 

IV 

input 

Vehicle  type:  2  for  booster,  4  for  RV 

COMMON  /HJRAJ/ 

input 

T,  RNG,  ALT,  N 

Routines  called: 

TERPL 

SUBROUTINE  TRPTTJ 


This  routine  computes  the  laminar  to  turbulence  transition  point 
for  the  drag  calculation  in  the  coupled  mode. 


COMMON  /DR I DR/ 

output 

XTR 

COMMON  /BC/ 

i  n  pu  t 

TCDEG,  TCRAD,  SINTC,  C0STC,  RN, 

RB 

COMMON  / FSPC/ 

input 

XM3 

C0W0N  /LPNC/ 

input 

XME2 ,  REFT2 

COMMON  /TRNSC/ 

output 

LAM I NR,  RNRTR 

COMMON  / PTMC/ 

input 

PI 

REXMIN 

data 

Minimum  Reynolds  number  allowed 
sition  (5x10s) 

for  tran- 

Routines  called: 

FNFME,  FNFRE 
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FUNCTI0N  V0LUME  (T,  HT,  HM,  RM,  RB) 


This  routine 

defines 

the  dust  cloud  volume  (kft3). 

T 

input 

Time  after  burst  (sec) 

HT 

input 

Dust  cloud  top  height  (kft) 

HM 

input 

Dust  cloud  middle  height  (kft) 

RM 

input 

Dust  cloud  middle  radius  (kft) 

RB 

input 

Dust  cloud  base  radius  (kft) 

C0MM0N  /THREAT/ 

i  nput 

I 

C0MM0N  /STUFF/ 

input 

TL0 

PI 

data 

ir  (3.1415927) 

FUNCTI0N  WFPK0D  (DUMMY) 

This  routine 

defines  peak  overdensity  at  shock  front  in  g/cm3. 

DUMMY 

input 

Dummy  argument 

C0MM0N  /WFRT/ 

input 

0PPK 

RH0Z 

data 

Density  of  air  at  sea  level  (1 .22E-3  g/cm3) 

Routines  called: 

AIR 

ill 


FUNCTI0N  WFPK0P  (R) 


This  routine  defines  peak  overpressure  at  the  shock  front  in  dynes/ 


1  R 

input 

Shock  front  radius  (cm) 

j  AC 

5 

data 

Constant  for  curve  fit 

! 

data 

Constant  for  curve  fit 

V 

j  ASTAR 

data 

Constant  for  curve  fit 

FUNCTI0N  WFPKV  (DUMMY) 

This  routine  defines  peak  air  particle  velocity  behind  the  shock 
front  in  cm/ sec. 


DUMMY 

input 

Dummy  argument 

C0MM0N  /WFRT/ 

input 

0PPK,  0DPK 

RH0Z 

data 

Density  or  air  at  sea  level  ( 1  . 22 E-3  g/cm3 
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SUBROUTINE  WSTUFF  (W,  SH0B) 


i 

I 


i 

t 

I 

i 

j 

I 

s 


This  routine  defines  constants  for  use  throughout  the  dust  cloud 
routines. 


w 

input 

Burst  yield  (Mt) 

SH0B 

input 

Scaled  height-of-burst  (ft/kt^3) 

C0MM0N  /CflNVR/ 

input 

CVL ,  CVM ,  CVT 

C0MM0N  /THREAT/ 

input 

I 

output 

W0LD,  SH0LD,  H0B ,  KDUST 

C0MM0N  /STUFF/ 

input 

TER 

output 

PCT,  ABO,  AFB ,  AHM,  AHT,  ARM,  BBO,  BFB , 

BHM, 

BHT,  BRM,  HMP ,  HMS ,  HMO ,  HTP ,  HTS ,  HTO , 

RBC , 

RB0,  RBS ,  RFS,  RMS,  R2M,  TBC ,  TB0,  TFS , 
TF1 ,  TL0 ,  TPH,  T2M 

TFW, 

C0MM0N  /MASSY/ 

output 

Al,  A2 ,  A5 ,  A6 ,  Bl,  B2,  B5,  B6,  Tl,  T2, 
T5,  WD1 ,  WD2 ,  WD3 ,  WD4 ,  WD5 ,  WD6 

T4 , 

AA1 

data  \ 

AA2 

data 

AA3 

data 

! 

>  Curve  fit  constants 

AA4 

data  | 

f 

AA5 

data  ' 

AA6 

data  t 

1 

THIRD 

data 

One  third 

EPS 

data 

Round  off  epsi Ion  (5x1 0"7) 

SH2 

data 

Maximum  SH0B  for  near  surface  burst  (cratering) 

mass  loading.  Above  this  air  burst  mass 
loading  used.  (20.0) 

Routines  called: 

TBREAK 

L 
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FUNCTI0N  XL I M  (X,  N,  I,  JM) 


This  function  finds  the  maximum  or  minimum  value  of  a  one  dimen 

sional  array. 

X  input  Array  of  values  to  be  searched 

N  input  Number  of  values  in  X  (Maximum  of  50) 

I  input  0peration  mode  flag:  I  >  0  for  maximum, 

I  <  0  for  minimum 

JM  output  Location  in  array  X  of  the  extreme  value 
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2.3  COWON  SLOCK  DESCRIPTION 

As  with  all  computer  codes  of  any  size  and  complexity,  NWEM  has  many 
common  blocks,  45  to  be  specific.  Besides  the  usual  function  of  transfer¬ 
ring  information  from  one  routine  to  another,  some  common  blocks  of  the 
NWEM  computer  code  are  the  connectors  with  which  the  executive  routine 
links  the  separate  environmental  and  exclusion  region  contour  generating 
routines  together,  while  others  serve  as  connectors  which  link  routines 
within  the  separate  environmental  and  contour  generating  models. 

To  provide  the  user  with  the  ability  to  better  use  the  NWEM  computer 
code,  detailed  descriptions  of  the  common  blocks  with  respect  to  their 
variables  are  provided  on  pages  following  in  this  section.  These  descrip¬ 
tions  present  a  brief  common  block  title  or  usage,  the  variables  in  the 
order  they  appear  in  the  common  block  and  their  dimensional  extent  if  an 
array,  the  variable's  definition  or  use  and  any  units  or  values  stored  in 
these  variables  where  appropriate.  Certain  common  blocks  will  have  their 
variables  assume  different  alphanumeric  names  in  different  subroutines. 

Only  one  name  is  presented  for  each  variable  in  the  common  block  definition; 
that  one  most  representative.  However,  all  common  blocks  of  this  type  will 
agree  for  variable  and  array  length  exactly.  No  overlay  of  variables  ex¬ 
tent  exists  for  different  routines. 
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C0MM0N  /BC/ 


Cone  Angle  and  Parameter  Storage  for  RV  Trajectory  and 
Drag  Calculation 


TCDEG 

RV  cone  angle  (deg) 

TCRAD 

RV  cone  angle  (radian) 

SINTC 

Sine  of  the  RV  cone  angle 

C0STC 

Cosine  of  the  RV  cone  angle 

TANTC 

Tangent  of  the  RV  cone  angle 

RN 

RV  nose  radius  (in) 

RB 

RV  base  radius  (in) 

$K 

RV  surface  roughness  (in) 

C0MM0N  /BOAT/ 

Blast  Footprint  Data 

TB(50) 

Time  after  burst  (sec) 

BL(50) 

Rectangular  footprint  half-length  (kft) 

BW(50) 

Rectangular  footprint  half-width  (kft) 

BC(50) 

Oownrange  position  (RV)  or  uprange  position  (booster) 
of  rectangular  footprint  center  (kft) 

BRMX(50) 

Annular  footprint  maximum  radius  from  burst  point  (kft) 

BRMN(50) 

Annular  footprint  minimum  radius  from  burst  point  (kft) 

PHI (50) 

Uprange  (RV)  or  downrange  (booster)  annular  segment 
internal  angle  (degrees) 

PH2(50) 

Downrange  (RV)  or  uprange  (booster)  annular  segment 
internal  angle  (degrees) 

NR 

Number  of  data  points  in  the  rectangular  footprint 
arrays  (Maximum  of  50) 

NC 

Number  of  data  points  in  the  annular  footprint 
arrays  (Maximum  of  50) 

no 


C0MM0N  /  BST/ 


Booster  Launch  Exclusion  Region  Data 


F0RADI  Static  region  radius,  centered  about  ground  zero  (NMi) 

F0TIMS  Time  duration  of  static  region  -  a  positive  quantity 

measuring  times  prior  to  burst  (sec) 

BlANG(5,7)  Blast  exclusion  contour  data  comprised  of  7  sequential 

time  sets  of  the  following  quantities  expressed  in 
AMM  nomenclature: 

1=1  -  BLATIM(7)  Blast  exclusion  region  time 

(negative  times  go  into  F0TIMS) 

(sec) 

1=2  -  BLA0UT ( 7 )  Outer  radius  of  the  blast  annular 
region  (NMi) 

1=3  -  SLAIN ( 7 )  Inner  radius  of  the  blast  annular 
region  (NMi) 

1=4  -  ANGNEA(7)  Total  annular  extent  of  the  near 
(uprange)  blast  sector  (deg) 

1=5  -  ANGFAR(7)  Total . annul ar  extent  of  the  far 
(downrange)  blast  sector  (deg) 

A  circular  region  is  expressed  by  setting  both  ANGNEA 
and  ANGFAR  equal  to  180. 

F0TLWC( 1,10)  Dynamic  exclusion  contour  data  (dust  and  ejecta/ pebble 

environments)  comprised  of  10  sequential  time  sets  of 
the  following  quantities,  expressed  in  AMM  nomenclature: 
1=1  -  F0TIME (10)  Dynamic  exclusion  region  time  -- 
only  positive  times  are  considered 
(NMi  ) 

1=2  -  F0LENT(1O)  Half-length  of  the  rectangular 
region  (NMi) 

1=3  -  F0WIDT(1O)  Half-width  of  the  rectangular 
region  (NMi) 

1=4  -  F0CENT(1O)  Distance  from  the  burst  ground  zero 
to  the  dynamic  region  center 
(positive  downrange)  (NMi) 

NBL  Number  of  blast  times  containing  non-zero  data  (maxi¬ 

mum  of  7).  Not  written  on  TAPE! 6 

NF0  Number  of  dynamic  region  times  containing  non-zero 

data  (maximum  of  10).  Not  written  on  TAPE16 

{ 


C0MM0N  /BURST/  Dust  Cloud  Parameters 


Dust  cloud  top  height  (kft) 

Dust  cloud  middle  height  (kft) 

Dust  cloud  middle  radius  (kft) 

Dust  cloud  base  radius  (kft) 

Dust  cloud  density  (g/cm3) 

Second  RV  burst  height  or  booster  launch  pt  height  (kft) 


C0MM0N  /C0NST/  Constants  for  Use  in  Blast  Routines 


C0MM0N  /C0NVR/ 


One  third 

Atmospheric  pressure  at  sea  level  (dynes/cm2) 
Atmospheric  sound  speed  at  sea  level  (cm/sec) 
Atmospheric  density  at  sea  level  (g/cm3) 
Atmospheric  temperature  at  sea  level  (°K) 


Unit  Conversion  Constants 


Length  scale  conversion  constant  (1.0  kft/kft) 

Mass  scale  conversion  constant  (1.0  kT/kT) 

Time  scale  conversion  constant  (60.0  sec/min) 

Density  scale  conversion  constant  (3.2036927x10  5 
g/cm3/kt-kft3) 


C0MM0N  / DDAT/ 


Oust  and  Pebble  Footprint  Data 


TD( 50 ) 

DL(50) 

DW(50) 

DC (50) 

NO 

C0MM0N  /DIMEN/ 

NY 

NX 

YMAX 

DY 

DX 

C0MM0N  /DRAG/ 
XCN(250) 

XCA(250) 

XA1 (10) 
XAM(25) 

W0A 

NA 

NM 


Time  after  burst  (sec) 

Rectangular  footprint  half-length  (kft) 

Rectangular  footprint  half-width  (kft) 

Downrange  position  (RV)  or  uprange  position  (booster) 
of  rectangular  footprint  center  (kft) 

Number  of  data  points  in  the  rectangular  footprint 
arrays  (Maximum  of  50) 

Dust  Erosion  Calculation  Grid  Parameters 

Number  of  uprange/downrange  positions  (100) 

Number  of  crossrange  positions  (23) 

Total  uprange/downrange  length  of  grid  (kft) 
Uprange/downrange  grid  interval  size  (kft) 

Crossrange  grid  interval  size  (kft) 

Drag  Data  Storage 

Normal  drag  force  coefficients  corresponding  to  Mach 
number  values  and  angle  of  attack  values 

Axial  drag  force  coefficients  corresponding  to  Mach 
number  values  and  angle  of  attack  values 

Angle  of  attack  values  (deg) 

Mach  number  values 

Weight  of  RV  divided  by  reference  area  of  RV  (lb/ft2) 
Number  of  angle  of  attack  values 
Number  of  Mach  number  values 


C0W0N  /DRIDR/ 
ITJF 

XTR 

ISTPAS 

C0MM0N  /DTH/ 
NAME (10) 


LABEL (10) 


C0MM0N  /EC0NS/ 


Drag  Calculational  Flag  Storage 

Trajectory  flag  indicating  that  this  is  a  coupled 
trajectory  drag  coefficient  calculation 

-Ratio  of  turbul ent/ Laminar  transition  point  along 
RV  body  to  the  R V  nose  radius 

Flag  to  initialize  variables 


Hollerith  Data  for  Program  Output  Label 


Environment  names 

1  -  NEUT  (Neutron  fluence) 

2  -  PGAM  (Prompt  gamma  dose  rate) 

3  -  XRAY  (X>ray  fluence) 

4  -  THML  (Thermal  fluence) 

5  -  EMP  (Peak  EMP  field  strength) 

6  -  BLST  (Peak  blast  loading  (either  net  or 

lateral  g's) 

7  -  DUST  (Dust 

8  -  EJCT  (Ejecta  and  pebbles) 

Subroutine  names  --  stored  as  two  4  character  words 
per  name,  left  justified 
1 ,2  EMPN 
3,4  STATIC 
5,6  STATRV 
7,8  RADI AT 
9,10  PT4 
11,12  PT1 7 
15,16  FENV 
17,18  MRV3 


Constants  for  Ejecta  Cloud  Dimensions 


BBB,  HIN,  HPK,  RIN,  RPK,  THP , 
THT,  TRP,  TRT 


Curve  fit  constants  for  ejecta  cloud 
dimensions  vs.  time 


C0MM0N  /EDAT/ 


Ejecta  Footprint  Data 


TE(50) 

EL ( 50 ) 

EW(50) 

EC ( 50 ) 

NE 

C0MM0N  /FRAT/ 
FM  ( 6 ) 

C0MM0N  / FSPC/ 

H2 

V2 

P8 

T8 

A8 

H8 

HT 

RH08 

XMU8 

XM8 

REFT8 

P8PT8 

Q8 


Time  after  burst  (sec) 

Rectangular  footprint  half-length  (kft) 

Rectangular  footprint  half-width  (kft) 

Downrange  position  (RV)  or  uprange  position  (booster) 
of  rectangular  footprint  center  (kft) 

Number  of  data  points  in  the  rectangular  footprint 
arrays  (Maximum  of  50) 

Blast  Kill  Criteria 

Blast  kill  criteria.  Locations  1-3  for  booster. 
Locations  4-6  for  RV 

Free  Stream  Parameter  Storage 

RV  altitude  at  beginning  of  time  step  (ft) 

-RV  velocity  at  beginning  of  time  step  (ft/sec) 

Free  stream  pressure  (lbf/ft2) 

Free  stream  temperature  (°R) 

Free  stream  sound  speed  (ft/sec) 

Free  stream  static  enthalpy  (BTU/lbm) 

Total  enthalpy  (BTU/lbm) 

Free  stream  density  (slug/ft3) 

Free  stream  viscosity  (slug/ft-sec) 

Mach  number 

Free  stream  Reynolds  number  per  foot  (ft-1) 

Free  stream  pressure  ratio  variables 
Free  stream  dynamic  pressure  (slug/ft-sec2) 
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C0MM0N  / GPC/ 


Atmosphere  and  Physical  Constants  for  Drag  Calculation 


CP8  Specific  heat  at  constant  pressure  (0.24  BTU/lbm  °R) 

G  Conversion  factor  from  slugs  to  Ibm  (gc)  (32.174) 

GAM8  Free  stream  ratio  of  specific  heats  (y)  (1.4) 

XJ  Conversion  factor  from  ft-lbf  to  BTUs  (778.0) 

PR  Prandfl  number  of  free  stream  air  (0.71) 

RGAS  Air  gas  constant  (R)  (1716  ft-lbf/sl ug-°R) 


C0MM0N  /HTRAJ/  Vehicle  Trajectory  in  HYDR0  Units 

T2 ( 75 )  Time  array:  Time  after  launch  for  booster,  Time 

before  impact  for  RV  (sec) 

R2 ( 75 )  Ground  range  from  launch  point  or  impact  point  (kft) 

A2(75)  Altitude  of  vehicle  (kft) 

V2 (75 )  Vehicle  velocity  (ft/sec) 

AN2 ( 75 )  Flight  path  angle  from  local  horizontal  (radians) 

N2  Number  of  data  points  in  trajectory  arrays 

(Maximum  of  75) 


C0MM0N  /HTRNS/  Hollerith  Data  Input  and  Passed  to  AMM 

LTYPE1  Identifier  for  Booster  or  incoming  RV 

LTYPE2  Identifier  for  detonating  RV 

LH0B  Height  of  burst  flag 


11.6 


C0MM0N  / ITRNS/  Integer  Input  Data 


ITYPE 

Flight  phase  flag  indicating  Booster  or  RV 

NTRJ 

Trajectory  input  flag 

NTRJX 

Absolute  value  of  trajectory  input  flag  | NTRJ | 

NDG1 

Number  of  input  Mach  numbers  and  associated  drag 
data  for  RV  trajectory  calculation 

NDG2 

Number  of  input  RV  angles  of  attack  and  associated 
drag  data  for  RV  trajectory  calculation 

IMRV 

MRV  flag 

INEUT 

Neutron  transmission  data  flag 

INET 

Absolute  value  of  neutron  transmission  data  flag, 

| INEUT) 

ixray 

X-ray  or  Gamma  ray  transmission  data  flag 

IXM 

Absolute  value  of  X-ray  or  Gamma  ray  transmission 
data  flag,  j IXRAY) 

C0MM0N  /LPNC/ 

Boundary  Layer  Property  Storage 

XMEW 

Local  Mach  number  at  the  outer  boundary  layer  edge 

REFTN 

Local  Reynolds  number  per  foot  at  the  outer  boundary 
layer  edge  (ft-1) 

RH0EW 

Local  density  at  the  outer  boundary  layer  edge  (slug/ ft3) 

TEW 

Local  temperature  at  outer  boundary  layer  edge  (°R) 

VEW 

Local  velocity  at  the  outer  boundary  layer  edge  (ft/sec) 

i 


C0MM0N  /MASSY/  Case  Constants  for  Dust  Cloud  Mass  Loading  Time  History 

t 


[ 


Al,  A2,  A3,  A4 ,  A5, 

62 ,  83,  84,  85,  B6, 

T2,  T3,  T4 ,  T5,  T6, 

WD3,  WD4,  WD5,  WD6 


A6 ,  Bl, 
TO,  Tl, 
WD1,  WD2 , 


Storage  arrays  for  dust  loading  curve 
fit  constants.  All  have  four  storage 
locations  for  four  burst  cases.  Each 
burst  case  identified  by  parameters 
in  C0MM0N  /THREAT/ 
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C0MM0N  /PUT/ 


Weapon  Output  and  Radiation  Transport  Data 


sx 

X-ray  source  strength  for  a  particular  yield  and  out¬ 
put  energy  fraction  (cal) 

TX(60) 

X-ray  fluence  buildup  (air  mass  transmission)  factor. 
Ratio  of  fluence  to  exoatmo spheric  fluence  for  a  par¬ 
ticular  air  mass  integral 

RTX ( 60 ) 

Air  mass  integral  values  corresponding  to  the  TX  array 
( gm/ cm2 ) 

$N 

Neutron  source  strength  for  a  particular  yield  and  out¬ 
put  number  (neutrons) 

TN(60) 

Neutron  fluence  buildup  factor.  Ratio  of  fluence  to 
exoatmospheric  fluence  for  a  particular  in  mass  integral 

RTN ( 60 ) 

Air  mass  integral  values  corresponding  to  the  TN  array 
( gm/ cm2 ) 

SG 

Prompt  gamma  source  strength,  expressed  as  a  dose  re¬ 
sponse  function,  for  a  particular  yield  and  output 
energy  fraction  (rad (Si J-cm2) 

TG(60) 

Prompt  gamma  dose  buildup  factor.  Assumed  identical 
to  the  energy  fluence  buildup  factor.  Ratio  of  energy 
fluence  to  exoatmospheri c  fluence 

RTG(60) 

Air  mass  integral  values  corresponding  to  the  TG  array 

ST 

Thermal  radiation  source  strength  for  a  particular 
yi el d  (cal ) 

N0X 

4 

Number  of  points  in  the  TX  and  RTX  arrays  containing 
data 

N0N 

Number  of  points  in  the  TN  and  RTN  arrays  containing 
data 

N0G 

Number  of  points  in  the  TG  and  RTG  arrays  containing 
data 

C0MM0N  /PTMC/ 


Constant  Storage  for  RV  Trajectory  Calculation 


PI 

XRAD 

XM0 

R0 

C0MM0N  /RAD/ 
CRIT (10) 

W0UT (4 ) 

ABURST 

PI 

AUPDN 

ALAST 

I SK I P ( 4 ) 

I PHASE 


ir  (3.141  592653589793) 

Radians  to  degrees  conversion  factor 

Gravitation  constant  times  the  mass  of  the  earth 
(1 .4076452x1 015  ft3sec2) 

Radius  of  earth  (2 .0925696x1 07  ft) 


Data  Transfer  for  Prompt  Radiation 

Criterion  used  to  define  the  lethal  volume  for  each 
environment.  These  are  input  data  paired  through  the 
SUBROUTINE  RADIAT  argument  list 
1=1  -  Neutron  fluence  (n/cm2) 

1=2  -  Peak  garrma  dose  rate  ( rad  (Si )/ sec  ) 

1=3  -  X-ray  fluence  (cal/cm2) 

1=4  -  Thermal  fluence  (cal/cm2) 

Weapon  source  strength  expressed  in  units  such  that 
division  by  an  area,  4Trr2,  corresponding  to  a  radius, 
r,  yields  a  quantity  which  can  be  directly  compared 
with  CRIT 

Burst  altitude,  an  input  quantity  passed  through  the 
SUBROUTINE  RADIAT  argument  list  (km) 

3.14159265 

Geometry  multiplier  for  finding  lethal  volume  dimen¬ 
sions  (+1.  for  a  calculation  above  the  burst  altitude, 
-1.  for  one  below  the  burst  altitude) 

Positive  angle  from  horizontal  to  the  maximum  radius 
ray  (normal  to  the  trajectory)  which  terminates  at 
the  lethal  volume  surface  point  being  found  (radians) 

Flag  for  each  environment,  set  equal  to  1  if  the 
environment  is  to  be  skipped  for  the  flight  phase 
being  addressed  or  if  a  criterion  is  not  input  (en¬ 
vironment  order  as  for  CRIT) 

Flight  phase  index  —  defined  by  data  passed  from 
input  through  the  SUBR0UTINE  RADIAT  argument  list 
2  -  Boost  Phase 
4  -  Reentry  Phase 


C0MM0N  /RAD/  (Con't) 


C0MM0N  /  RR1  / 


T(7 ,4) 


RH(7,4) 


R(7,4) 


RDEL( 7 ,4 ) 


RUP(7,4] 


RDN(7,4) 


A( 7,4 ) 
PERCNT(7,4) 


Index  for  determining  a  reasonable  number  of  itera¬ 
tions  to  begin  converging  on  the  lethal  volume  size 
in  subroutines  PT4  and  PT1 7 

Index,  I,  of  the  time  array  (T(I),  1=1,7)  correspond 
ing  to  the  maximum  range  as  calculated  in  SUBR0UTINE 
PTMAX  so  that  the  proper  time  index  for  the  final 
point  calculated  by  SUBR0UTINE  PTLAST  can  be  found 


Parameters  Defining  the  Radiation  Environment  Lethal 
Vo  1  ume 

Time  of  flight  between  lethal  volume  point  (i.e.,  a 
hardness  plane  in  the  lethal  volume)  and  the  earth's 
surface  (sec) 

Horizontal  lethal  volume  radius  at  time,  T.  This  is 
equal  to  the  instantaneous  exclusion  region  radius 
on  the  ground  (km) 

Radius  from  burst  to  lethal  volume  surface  correspond¬ 
ing  to  time,  T  (km) 

Horizontal  displacement  of  trajectory  between  lethal 
volume  point  corresponding  to  T  and  the  surface.  This 
is  positive  in  the  downrange  direction  (km) 

Maximum  uprange  extent  of  the  surface  exclusion  region 
corresponding  to  T,  a  negative  quantity.  Measured 
from  burst  ground  zero  (km) 

Maximum  downrange  extent  of  the  surface  exclusion 
region  corresponding  to  T,  a  positive  quantity.  Meas¬ 
ured  from  burst  ground  zero  (km) 

Altitude  corresponding  to  Time,  T  (km) 

Diagnostic  information  —  the  percentage  error  in 
convergence  of  the  environment  to  the  desired  environ¬ 
ment  criterion,  CRIT  at  the  lethal  volume  surface  point 
corresponding  to  T 


(a)Each  parame* ar  in  this  common  block  is  dimensioned  up  to  ITM(I)  =  ITMAX, 
a  maximum  or  7  times  and  NENV,  a  maximum  of  4  environments.  Time  index 
1  corresponds  to  the  top  of  the  lethal  volume. 


C0MM0N  / RR2/ 
TEX(7)^a) 


RHEX( 7 ) 
REX(7) 
RDELEX(7) 
RUPEX ( 7 ) 

AX(7) 
RDNEX(7 ) 

TMIN 

TMAX 

RHMAX 

RUPMIN 

RDNMAX 

NAMEX(7) 

NAMET1 

NAMET2 


Extreme  Values  of  Radiation  Environment  Lethal  Volume 
and  Exclusion  Region  Parameters 

A  set  of  exclusion  region  times  defined  to  span  the 
total  range  of  times,  T,  for  all  radiation  environ¬ 
ment.  Converted  to  negative  times  for  the  boost 
phase  (sec) 

Maximum  horizontal  radius,  RH,  for  all  radiation  en¬ 
vironments  (km) 

Range,  R,  for  environment  maximizing  RH  at  each 
time  (km) 

Horizontal  displacement,  RuEL,  for  environment  maxi¬ 
mizing  RH  at  each  time  (km) 

Uprange  contour  extent,  RUP,  for  environment  maximiz¬ 
ing  RH  at  each  time.  Measured  from  burst  ground 
zero  (km) 

Altitude  corresponding  to  TEX  (km) 

Downrange  contour  extent,  RDN,  for  environment  maxi¬ 
mizing  RH  at  each  time.  Measured  from  burst  ground 
zero  (km) 

Minimum  value  of  TEX ( I )  for  the  start  of  the  exclusion 
region  (the  most  negative  time  in  the  boost  phase) 
(sec) 

Maximum  value  of  TEX(I)  for  the  end  of  the  exclusion 
region  (sec) 

Maximum  value  of  RHEX(I) 

Minimum  value  of  RUPEX(I) 

Maximum  value  of  RDNEX(I) 

Name  of  radiation  environment  maximizing  RH  at  each 
time  --  drawn  from  NAME(I)  (hollerith  variable) 

Radiation  environment  name  corresponding  to  TMIN 
(hollerith  variable) 

Radiation  environment  name  corresponding  to  TMAX 
(hollerith  variable) 


NAMEH 


Radiation  environment  name  corresponding  to  RHMAX 
(hoi lerith  variable) 
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C0MM0N  /RR2/  (Con't) 


NAMEU  Radiation  environment  name  corresponding  to  RUPMIN 

(hoi lerith  variable) 

NAMED  Radiation  environment  name  corresponding  to  RDNMAX 

(hollerith  variable) 

(a)The  number  of  points  calculated  is  MAXITM  (See  C0MM0N  / RRR/ ) .  The 
first  time  index  corresponds  to  the  top  of  the  lethal  volume. 


C0MM0N  /RRR/  Radiation  Exclusion  Region  Parameters 


UMAX 

TIMIN 

JINDEX 

MAXITM 

IL0OPS 

ITM(4) 


Maximum  time  interval,  lethal  volume  to  surface  for 
all  environments,  a  positive  quantity  (sec) 

Minimum  time  interval,  a  lethal  volume  to  surface  for 
all  environments,  a  positive  quantity  (sec) 

Index  of  the  first  radiation  environment  having  valid 
lethal  volume  data.  Ranges  from  1  to  4  wi th  0  indica¬ 
ting  no  valid  data  (environment  types  are  as  defined 
by  NAME(I),  1=1,4  in  C0MM0N  / DTH/ ) 

Maximum  number  of  times  specified  for  any  of  the  en¬ 
vironments.  Ranges  from  4  to  7  if  any  data  are  valid. 
If  two  or  more  radiation  environments  are  calculated, 
this  redefined  as  7  in  SUBR0UTINE  RADMAX 

The  running  number  of  environments  for  which  lethal 
volumes  have  been  sized.  Ranges  from  0  to  4 

Number  of  times  calculated  for  each  environment. 

Ranges  from  4  to  7 
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C0MM0N  / R VC/  Reentry  Vehicle  Targeting  Exclusion  Region  Oata 

FRSLEN  Half-length  of  the  elliptical  static  region  (NMi) 

FRSWID  Half-width  of  the  elliptical  static  region  (NMi) 

FRSCEN  Distance  from  the  burst  ground  zero  to  the  static 

region  center  (NMi) 

TIMST  Duration  of  the  static  region  (sec) 

TMFLWC (1,10)  Dynamic  exclusion  contour  data  (blast,  dust,  and 

ejecta/pebble  environments)  comprised  of  10  sequen¬ 
tial  time  sets  of  the  following  quantities,  expressed 
in  AMM  nomenclature: 

1=1  -  TIMDY(IO)  Dynamic  exclusion  region  time  (min) 
1=2  -  FRLENT(IO)  Half-length  of  the  rectangular 
region  (NMi) 

1=3  -  FRWIDT(IO)  Half-width  of  the  rectangular 
region  (NMi) 

1=4  -  FRCENT(IO)  Distance  from  the  burst  ground 

zero  to  the  dynamic  region  center 
(NMi) 

NTM  Number  of  dynamic  region  times  containing  non-zero  data 

(maximum  of  10).  Not  written  on  TAPE! 6 


C0MM0N  / S0R/  Stored  Radiation  Transmission  Data 

XX (2 )  Value  of  X-ray  energy  fraction 

TX I (60 , 2 )  X-ray  energy  fluence  build  up  factor  normalized  to 

exoatmospheri c  fluence 

RTX 1(60,2)  Air  mass  integral  corresponding  to  X-ray  energy 

fluence  build  up  factors  (gm/cm2) 

XN(1)  Values  of  neutron  output 

TN I ( 60 , 1 )  Neutron  number  of  fluence  build  up  factor  normalized 

to  exoatmospheric  fluence 

RTN I ( 60 , 1 )  Air  mass  integral  corresponding  to  neutron  number 

fluence  build  up  factor  (gm/cm2) 

XG(1)  Value  of  gamma  ray  energy  to  dose  conversion  factor 

TGI ( 60 , 1 )  Prompt  ganma  ray  energy  build  up  factor  normalized 

to  the  exoatmospheric  fluence 


C0MM0N  /S0R/  (Con ' t) 


RTGI(60,1 ) 


NQX(2) 


NQN(1 ) 


NGN(1 ) 


Air  mass  integral  build  up  factor  corresponding  to 
prompt  gamma  ray  energy  build  up  factor  (gm/cm2) 

Number  of  data  points  for  each  stored  X-ray  trans¬ 
mission  data  set 

Number  of  data  points  for  each  stored  neutron  trans¬ 
mission  data  set 

Number  of  data  points  for  each  stored  gamma  ray 
transmission  data  set 

Number  of  stored  X-ray  transmission  data  sets 
Number  of  stored  neutron  transmission  data  sets 
Number  of  stored  ganrna  ray  transmission  data  sets 


C0MM0N  /SPC/ 


Static  Pressure  Storage  for  Drag  Calculation 


C0MM0N  /STA/ 


TMN(2) 


TMX(2) 


RHMX(2) 


RUPMIN(2) 


Static  pressure  at  boundary  layer  edge  (lb/ft2) 
Static  pressure  to  free  stream  pressure  ratio 


Extreme  Values  of  the  Blast  and  Radiation  Exclusion 
Region  Contours  Over  All  Times 

Static  exclusion  contour  start  time;  negative  for 
boost,  positive  for  reentry  (sec) 

Static  exclusion  contour  ending  time;  zero  for  boost, 
positive  for  reentry  (sec) 

Maximum  instantaneous  radius  of  the  surface  exclusion 
region.  This  is  also  the  maximum  cross  range  contour 
radius  (km) 

Maximum  exclusion  contour  uprange  extent  (more  nega¬ 
tive  direction)  (km) 


(a)The  first  index  corresponds  to  radiation  data,  the  second  to  blast 
(since  only  blast  extends  to  significant  negative  times  with  radiation 
contours  for  the  boost  phase).  All  names  with  index  equal  to  2  are 
"BLST"  in  SUBROUTINE  STATIC  for  the  boost  phase  and  are  not  used  by 
SUBROUTINE  STATRV  for  the  reentry  phase. 


C0MM0N  /STA/  (Co n't) 

RDNMX(2)  Maximum  exclusion  contour  downrange  extent  (more 

positive  direction)  (km) 


NMT1(2)  Name  of  environment  corresponding  to  TMN  (hollerith 

variable) 

NMT2(2)  Name  of  environment  corresponding  to  TMX  (hollerith 

variable) 

NMRH(2)  Name  of  environment  corresponding  to  RHMX  (hollerith 

variable) 

NMRU(2)  Name  of  environment  corresponding  to  RUPMN  (hollerith 

variable) 

NMRD(2)  Name  of  environment  corresponding  to  RDNMAX  (hollerith) 

variable) 

IDATA(2)  Flag  indicating  presence  of  valid  exclusion  contour 

data  for  environment  type  I  of  IDATA  is  .NE  -  0 
1=1  -  Prompt  radiation 
1=2  -  Blast 


C0MM0N  /STUFF/  Case  Constants  for  Dust  Cloud  Dimension  Time  Histories 


PCT , 

ABO, 

AFB , 

AHM, 

AHT , 

ARF 

ARM, 

BBO, 

BFB , 

BHM, 

BHT , 

BRF 

BRM, 

HMP , 

HMS , 

HMO , 

HTP , 

HTS 

HTO, 

RBC, 

RB0, 

RBS , 

RFS , 

RF1 

RMS, 

R2M, 

TBC, 

TB0, 

TBS, 

TER 

TFS, 

TFW, 

TF1, 

TL0, 

TPH, 

TRS 

T2M 

RFB 


Storage  arrays  for  cloud  dimension 
curve  fit  constants.  All  have  four 
storage  locations  for  four  burst  cases. 
Each  burst  case  identified  by  para¬ 
meters  in  C0MM0N  /THREAT/ 


Fireball  radius  as  defined  in  EARLY 
(kft) 


C0MM0N  /THER/  Default  Values  for  Input  Quantities 
THER1 


DUSTFL 


Thermal  transmission  default  value 
Dust  cut-off  time  default  value  (25  min) 


C0MM0N  /THREAT/  Identifiers  for  Burst  Conditions  Which  A1 ready  Have 
Constants  Stored 


W0LD(4) 

SH0LD(4) 

H0B(4) 

KDUST(4) 


Index  of  current  case 

Number  of  conditions  stored  (Maximum  of  4) 

Burst  yield  array  (Mt) 

Scaled  height  of  burst  array  (ft/kt1,/3) 

Height  of  burst  array  (kft) 

Dust  cloud  flag: 

0  -  no  dust  cloud 
1  -  normal  dust  cloud 


C0MM0N  /TRAJ/ 


TI  (75) 


R1  (75) 
ALT (75) 
VI (75) 
A1  (75) 


NTRAJ 


QT2 (75,1 ) 
QR2(75,1) 
QAL2 (75,1) 
QV2 (75 ,1 ) 
QA2 (75,1) 
NT  ( 1  ) 


Vehicle  Trajectory  Data  Storage 
Time  array 

Booster  -  time  after  launch  (sec) 

RV  -  time  before  impact  (sec) 

Ground  range  from  launch  point  or  to  impact  point  (NMi) 

Altitude  (ft) 

Velocity  (ft/ sec) 

Flight  path  angle  with  respect  to  local  horizontal 
(deg) 

Number  of  data  points  in  trajectory  arrays  (75  maximum) 


C0MM0N  /TRAJX/  Stored  Booster  Trajectory  Data 


Trajectory  time  data  (sec) 

Trajectory  range  data  (nmi) 

Trajectory  altitude  data  (ft) 

Trajectory  velocity  data  (ft/sec) 

Trajectory  booster  angle  data  (deg) 

Number  of  data  points  for  each  stored  trajectory 
Number  of  stored  trajectories 


C0MM0N  /TRNS/ 


Input  Data  Storage 


BLST1 

BLST2 

BLST3 

DUST! 

DUST2 

DUST3 

PEBB1 

PEBB2 

THER1 

THERF 


Blast  vulnerability  criteria  data 

Booster  -  free  field  overpressure  (psi) 

RV  -  total  acceleration  (g's) 

Blast  vulnerability  criteria  data 

Booster  -  ratio  of  overpressure  to  ambient  pressure 
RV  -  axial  aceleration  (g's) 

Blast  vulnerabil i ty  criteria  data 

Booster  -  dynamic  pressure  times  angle  of  attack 
(psi -deg) 

RV  -  normal  acceleration  (g's) 

Dust  vulnerability  criteria,  intercepted  kinetic 
energy  (kJ/cm2) 

Dust  vulnerability  criteria,  cloud  cutoff  time  (min) 

Cloud  cutoff  time  used  for  calculation  (min) 

Pebble/ejecta  vulnerability  criteria,  critical  par¬ 
ticle  diameter  (cm) 

Pebble/ejecta  vulnerability  criteria,  critical  hit 
density  (#/cm2) 

Thermal  vulnerability  criteria,  thermal  fluence 
(cal/crn2) 

Fluence  multiplication  factor  for  transmission  cor¬ 
rection 


XNEUT1 

XGM1 

EMP1 

RVHT 

RVANG 

RVVEL 

RVMAS 

RVARA 


Neutron  vulnerability  criteria,  neutron  fluence 
(n/cm2 ) 

X-ray/gaiuna  ray  vulnerability  criteria 
Booster  -  X-ray  energy  fluence  (cal /cm2 
RV  -  Peak  gamma  ray  dose  rate  (rad(Si )/sec) 

EMP  vulnerability  criteria,  peak  EMP  field  strength 
(V/m) 

RV  initial  reentry  altitude  (ft) 

RV  initial  reentry  angle  (deg) 

RV  initial  reentry  velocity  (ft/sec) 

RV  mass  (slugs) 

RV  reference  (base)  area  (ft2) 
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C0MM0N  /TRNS/  (Con't) 


RVNR 

RVBR 

RVCHA 

RVSR 

RVCYHA 

RVFHA 

RVCL 

RVCYC 

RVFL 

YLD 

YH0B2 

XNEUTX 

XFRAC 

GAMDS 

GAMPL 

RADMAX 

RADTIM 

XTM1 

XTM2 

TIMX(75) 
ALT(75) 
RNG( 75) 
VEL (75 } 
ANG(75) 


RV  nose  radius  (in) 

RV  base  radius  (in) 

RV  cone  half  angle  (deg) 

RV  surface  roughness  (in) 

RV  cylinder  half  angle  (deg) 

RV  flare  half  angle  (deg) 

RV  cone  length  along  axis  (in) 

RV  cylinder  length  along  axis  (in) 

RF  flare  length  along  axis  (in) 

Warhead  yield  of  detonating  RV  (kt) 

Burst  altitude  of  detonating  RV 

Neutron  output  (n/kt) 

X-ray/gamma  ray  energy  fraction 

Ganna  ray  energy  to  dose  conversion  factor 
(rad(Si )cm2/cal ) 

Gamma  ray  pulse  width  (nsec) 

MRV  pattern  radius  (NMi) 

Time  between  bursts  (sec) 

Ganna  ray  vulnerability  criteria  used,  peak  gaima 
ray  dose  rate  (rad(Si )/sec) 

X-ray  vulnerability  criteria  used,  x-ray  energy 
fluence  (cal/cm2) 

Booster  trajectory  data,  time  after  launch  (sec) 

Booster  trajectory  data,  altitude  (ft) 

Booster  trajectory  data,  ground  range  (NMi) 

Booster  trajectory  data,  velocity  (ft/ sec) 

Booster  trajectory  data,  flight  path  angle  above 
the  horizon  (deg) 


C0MM0N  /TRNS/ 

AMX(25 ) 

ANGX(IO) 

CN(250) 

CA(250) 

RH0RN (60 ) 

XNETX(60) 

RH0RX(6O) 

XPHI (60 ) 

zmm  /trnsc/ 

LAMINR 

RNRTR 

XKTR 

C0MM0N  /TRNSL/ 
NDG1M 

NDG2M 

NTRJXM 

INETM 


Ion ’t) 


Mach  number  values  for  drag  data 

Angle  of  attack  values  for  drag  data  (deg) 

Normal  drag  force  coefficients  corresponding  to 
the  Mach  number  values  and  angle  of  attack  values 

Axial  drag  force  coefficient  corresponding  to  the 
Mach  number  values  and  angle  of  attack  values 

Air  mass  integral  values  for  neutron  transmission 
data  (g/cm2) 

Normalized  neutron  transmission  values  as  a  function 
of  air  mass  integral 

Air  mass  integral  values  for  X-ray/ gamma  ray  trans¬ 
mission  data  (g/cm2) 

Normalized  X-ray/gamma  ray  transmission  values  as  a 
function  of  air  mass  integral 


Drag  Coefficient  Parameter  Storage 

Laminar  flow  flag  (0/1  -  laminar/ turbulent  or  mixed) 

Ratio  of  the  RV  nose  radius  to  the  laminar  turbulence 
transition  point 

Transition  point  parameter 


Storage  Limits  for  Input  Quantities 

Storage  limit  for  Mach  numbers  and  associated  drag 
data  for  RV  trajectory 

Storage  limit  for  RV  angles  of  attack  and  associated 
drag  data  for  RV  trajectory 

Storage  limit  for  booster  trajectory  data 

Storage  limit  for  neutron  transmission  data 

Storage  limit  for  X-ray  or  gamma  ray  transmission  data 


IXTRM 


C0MM0N  /VOLLY/ 


V 

C0MM0N  /WFRT/ 

PRAD 

0PPK 

0DPK 

VPK 

0PR 

0DR 

VR 

RZP 

RZD 

RZV 

0PMN 

0DMN 

VMN 


Oust  cloud  volume  (kft3) 

Air  Blast  Waveform  Parameters  for  1  Kt  Free  Air  Burst 
at  Sea  Level 

Shock  radius  (cm) 

Peak  overpressure  at  PRAO  (dynes/cn2) 

Peak  overdensity  at  PRAD  (g/cm3) 

Peak  air  velocity  at  PRAD  (cm/sec) 

Overpressure  (dynes/cm2) 

Overdensity  (g/cm3) 

Air  velocity  (cm/ sec) 

Range  to  end  of  overpressure  positive  phase  (cm) 
Range  to  end  of  overdensity  positive  phase  (cm) 
Range  to  end  of  air  velocity  positive  phase  (cm) 
Minimum  overpressure  (dynes/cm2) 

Minimum  overdensity  (g/cm3) 

Minimum  air  velocity  (cm/sec) 
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C0MM0N  / WPN/ 


Weapon  Radiation  Data 


W 

FN 

FG 

FX 

FT 

DELTG 

TFACT 

COMMON  / XTITL/ 
IEV0UT 

TITLE (13) 


Burst  yield  (kt) 

Neutron  output  source  strength;  set  equal  to  2.41  E+23 
if  not  input  (neutrons/kt) 

Gamma  output  energy  fraction;  set  equal  to  0.003  in 
the  unclassified  version  of  NWEM  if  not  input.  See 
Section  3.4  of  Volume  I  for  a  discussion  of  the  pro¬ 
per  value  for  specific  warheads 

X-ray  output  energy  fraction;  set  equal  to  0.75  if 
not  input 

Thermal  output  energy  fraction  for  low  altitude 
bursts.  Defined  by  a  yield-dependent  functional  re¬ 
lationship  if  not  input 

Effective  prompt  gamma  fall-width,  half-maximum 
(FWHM)  pulsewidth.  Set  equal  to  20.  if  not  input 
See  section  3.5  of  Volume  I  for  a  discussion  of  the 
proper  value  for  specific  warheads  (nsec).  Conver¬ 
sion  factor  from  gamma  energy  fluence  to  silicon 
dose,  valid  at  the  anticipated  range  (air  mass 
integral)  for  the  prompt  gamma  dose  rate  lethal 
volume  surface;  set  equal  to  1.4  E+03  if  not  input 
(rad(Si )/(cal/cm2) 

Factor  by  which  the  thermal  radiation  environment 
is  multiplied  to  account  for  atmospheric  transport 
and  cloud/earth  surface  albedo 


Output  Flag  and  Problem  Title  Storage 

Flag  indicating  level  of  output 

0  -  File  16  AMM  exclusion  contour  region  data 

1  -  In  addition  to  above,  input  and  AMM  exclusion 

contour  region  data  output 

2  -  In  addition  to  above  an  expansion  of  input 

data  results  and  exclusion  contour  region 
data  by  environment 

3  -  In  addition  to  above,  program  diagnostic/debug 

data  output 

Storage  for  problem  title 


2.4  BLOCK  DATA  DESCRIPTION 

The  block  data  that  is  in  the  present  version  of  the  NWEM  computer 
code  is  composed  of  tvo  types  of  variables.  This  division  of  the  vari¬ 
ables  into  two  groups  results  from  their  use  or  function.  The  two  uses 
of  the  variables  are  for  environment  model  and  program  functional  data 
and  for  prestored  input  related  data.  Those  variables  defined  as  en¬ 
vironment  model  variables  are  those  variables  which  are  associated  with 
the  model  descriptions  which  predict  the  environments.  Those  variables 
defined  as  program  functional  variables  are  those  variables  associated 
with  the  program,  such  as  storage  limitations  and  output  labels.  Those 
variables  defined  as  prestored  input  related  variables  are  used  to  store 
data  that  would  otherwise  be  read  as  input  from  the  input  file  (5).  The 
necessity  for  those  data  used  for  the  environment  models  and  program 
functional  data  is  evident.  Those  data  are  contained  in  the  commons 
/TRNSL/,  /THER/,  /DTH/,  /C0NVR/ ,  /THREAT/,  /MASSY/,  and  /STUFF/  and  are 
defined  in  their  respective  common  definitions  in  the  previous  section  on 
common  blocks. 

The  use  of  block  data  for  prestored  input  related  data  is  not  as 
evident.  At  the  time  the  NWEM  computer  code  was  being  designed,  it  was 
thought  that  for  code  check-out  the  most  straight  forward  method  of  hand¬ 
ling  input  related  data  would  be  using  block  data.  By  prestoring  that  part 
of  the  input  data  such  as  radiation  transmission  data  and  booster  trajec¬ 
tory,  the  amount  of  input  required  was  reduced  and  problem  input  facil¬ 
itated,  making  check-out  of  the  code  easier.  To  a  limit,  this  technique  of 
prestored  input  related  data  can  be  extended  to  include  more  and  more 
data  sets;  however,  a  limit  does  exist  (See  section  3.1). 

The  prestored  input  related  data  consists  of  two  types,  radiation 
transmission  data  and  booster  trajectory  data.  In  the  NWEM  computer  code, 
the  nuclear  radiation  environments  are  modeled  using  "build-up  factors", 
air  transmission  factors  which  are  the  value  of  the  air  transported  en¬ 
vironment  normalized  to  the  space  environment  value.  (See  section  3.4  of 
Volume  I).  The  ANISN  discrete  ordinates  code-  has  been  used  to  determine 
air  transmission  functions  for  representative  neutron,  X-ray  and  peak 
prompt  garrma  doses.  These  data;  a  representative  neutron  transmission 
set,  two  representative  X-ray  transmission  sets  and  one  representative 
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uncoil ided  prompt  gamma  set;  are  stored  in  the  variables  of  common  / S0R/ . 
(See  previous  section).  Since  booster  trajectories  are  not  calculated  in 
NWEM  and  are  required,  they  must  be  supplied.  One  representative  trajec¬ 
tory  is  stored  in  coimion  /TRAJX/  (See  previous  section). 


rvjg.,2svm 


2.5  INPUT  GUIDE 


Input  to  the  NWEM  computer  code  is  read  from  the  normal  input  file 
(5)  and,  at  present,  no  other  input  files  are  required.  The  input  data 
nas  been  grouped  into  two  major  sets,  those  parameters  associated  with  the 
vehicle  (booster  or  RV)  characteristi cs  necessary  to  define  the  exclusion 
region  criteria  and  those  parameters  associated  with  the  detonating  RV  to 
define  the  environments  resulting  from  the  nuclear  explosion.  Within  each 
major  group,  the  input  has  also  been  grouped  to  keep  parameters  of  like 
kind  together. 

For  the  booster  flyout  exclusion  region  contours,  three  sets  of  data 
are  needed.  The  first  data  set  is  the  booster  and  problem  identification 
and  flags.  While  booster  identification  is  not  needea  by  NWEM  for  calcula- 
tional  purposes  it  is  required  for  output  to  the  AMM  output  file  (16).  The 
second  set  of  data  is  the  vulnerability  criteria  data.  For  the  booster 
these  criteria  are  blast,  dust,  pebble/ejecta,  thermal,  neutron,  x-ray  and 
high  altitude  EMP.  The  third  set  of  data  is  the  booster  trajectory  data. 

In  the  present  version  of  th?  code  the  booster  trajectory  can  also  be 
specified  as  being  one  prestored  in  the  block  data. 

For  the  RV  fratricide  exclusion  region  contours,  again  three  sets  of 
data  are  needed.  The  first  data  is  the  RV  and  problem  identification  and 
flags.  As  with  the  booster,  the  RV  identification  is  for  output  to  the  AMM 
output  file  (16).  The  second  set  is  the  vulnerability  criteria  data  which 
for  RVs  are  blast,  dust,  pebble/ejecta,  neutron  and  gamma  ray.  The  third 
set  of  data  is  the  vehicle  parameters  to  compute  the  RV  trajectory.  This 
set  of  data  can  take  two  forms.  One  form  is  the  input  of  drag  force  coef¬ 
ficients  arrays  to  be  used  for  both  the  trajectory  calculation  and  the  RV 
blast  vulnerability  calculation.  The  other  form  is  the  input  of  additional 
RV  geometry  to  be  used  to  calculate  the  drag  force  coefficients  array  to 
be  used  for  both  the  trajectory  calculation  and  the  RV  blast  vulnerability 
calcul ation . 

The  parameters  associated  with  the  detonating  RV  are  grouped  into  two 
sets,  the  RV  identifiers  and  flags  and  the  weapon  characteristics  including 
the  air  transport  of  the  radiation.  Again,  there  are  some  RV  identifiers 
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which  are  only  needed  for  output  to  the  AMM  output  file.  And  in  the  pre¬ 
sent  version  the  air  transports  for  the  radiation  environments  can  be 
specified  as  being  prestored  in  block  data. 

Because  AMM  treats  MRVs  in  a  coupled  fashion  (see  Volume  I),  there 
are  data  required  as  part  of  the  RV  identifiers  specifying  MRV  type. 

In  the  case  of  the  MRV  type  in  which  the  RVs  are  in-line  and  having  dif¬ 
ferent  weapon  yields,  characteristics  or  burst  altitudes,  the  weapon  charac¬ 
teristics  cards  have  to  be  input  for  a  second  pass  of  the  code  to  get  the 
second  exclusion  region  contours.  (At  present  this  option  of  MRV  is  not 
impl emented) . 

Besides  grouping  the  data  consistently,  the  field  structure  of  the 
cards  has  been  formatted  to  facilitate  user  handling  of  the  input.  All 
formats  are  arranged  for  fields  to  be  5  or  10  characters  only.  No  other 
increments  are  used. 

On  the  following  pages  of  this  section  the  input  to  the  code  is  given 
and  explained  for  each  card  in  detail.  Restrictions,  defaults  and  other 
requirements  have  been  included  to  make  this  input  self  explanatory.  Ad¬ 
ditionally,  the  user  is  directed  to  the  sample  problems  for  further  ex- 
pi  anation. 


Card  1  Problem  Identification  and  Flags 

Format  15,  A4,  IX,  15,  4X,  AT,  13A4 

ITYPE  Flight  phase  flag 

Booster  -  2 
RV  -  4 

LTYPE1  Vehicle  identifier 


NTRJ  Trajectory  input  flag 

Booster  -  NTRJ  >,  0,  use  prestored  trajectory  set  number 
NTRJ  >  0,  input/NTRJ/  trajectory  data 
card  sets  as  Card  16 
RV  -  NTRJ  2  0,  calculate  drag  data 

-  NTRJ  <  0,  input  drag  data  as  Cards  10  through  14 

LEV0UT  Output  option  flag 

0  -  Only  AMM  exclusion  region  data  are  output  on  AMM 
output  file  (16) 

1  -  In  addition  to  above,  input  is  echoed  and  AMM  ex¬ 

clusion  region  data  are  output  to  normal  output 
file  (6) 

2  -  In  addition  to  above,  the  input  echo  is  expanded 

and  exclusion  regions  by  enviro,,.  ant  are  output  on 
normal  output  file 

3  -  In  addition  to  above,  more  detailed  exclusion  region 

data  are  output  on  normal  output  file 
Error  and  program  messages  are  always  output  on  normal 
output  file 

TITLE  Problem  title 


Card  2 
Format 


Blast  Vulnerability  Criteria  Data 
3E10.0 


BLST1  Booster  -  Free  field  overpressure  (psi) 

RV  -  Total  acceleration  (g's) 

BLST2  Booster  -  Ratio  of  overpressure  to  ambient  pressure 
RV  -  Axial  acceleration  (g's) 

BLST3  Booster  -  Dynamic  pressure  times  angle  attack  (psf-deg) 
RV  -  Normal  acceleration  (g's) 

If  all  3  values  are  0.,  no  blast  exclusion  region  contour  will  be 

calculated.  A  composite  contour  is  calculated  for  combinations  of 
the  cri teria . 


Card  3  Dust  Vulnerability  Criteria 
Format  2E10.0 

DUST1  Intercepted  kinetic  energy  (kj/cm2) 

DUST2  Cloud  cut-off  time  (min) 

If  both  values  are  0,  no  dust  exclusion  region  contour  will  be  cal¬ 
culated.  If  DUST1=0.  and  DUST2^0.,  cloud  avoidance  will  be  used.  For 
DUST2  <  0.  or  for  DUST2=0.  and  DUST!  >0.,  the  default  (25  min)  will 
be  used  for  cloud  cut-off  time. 
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Format 


Pebble/Ejecta  Vulnerability  Criteria 
2E10.0 

PEBB1  Critical  particle  diameter  (cm) 

PEBB2  Critical  hit  density  (#/cm2) 

If  both  values  are  0.,  no  pebble/ejecta  exclusion  region  contours  are 
calculated.  If  either  is  not  0.,  the  code  presently  only  does  cloud 
avoidance  to  a  predetermined  cut-off  time. 

Card  5  Thermal  Vulnerability  Criteria 

(Enter  only  if  ITypE  =  2,  Booster) 

Format  2E10.0 

THER1  Thermal  fluence  (cal /cm2) 

THERF  Fluence  multiplication  factor  for  transmission  correc¬ 

tions  (defaults  to  1.0) 

If  THER1=0.,  no  thermal  exclusion  region  contour  will  be  calculated. 

Card  6  Neutron  Vulnerability  Criteria 

Format  El  0.0 

XNEUT1  Neutron  Fluence  (n/cm2) 

If  XNEUT1 =0 . ,  no  neutron  exclusion  region  contour  will  be  calculated. 

Card  7  X-ray/Gamma  Ray  Vulnerability  Criteria 

Format  E10.0 

XMG1  Booster  -  X-ray  fluence  (cal/cm2) 

RV  -  Peak  gamma  ray  dose  rate  (rad(Si )/sec ) 

If  XMG1=0.,  no  X-ray/gamma  ray  exclusion  region  contour  will  be  cal¬ 
culated. 

Card  8  EMP  Vulnerabil ity  Criteria 

(Enter  only  if  ITYPE=2,  booster) 

Format  El  0.0 

EMP1  Peak  EMP  field  strength  (V/m) 

If  EMP1=0.,  no  exclusion  region  contour  will  be  calculated 

Card  9  Reentry  Vehicle  Parameters 

(Enter  only  if  ITYPE=4,  RV) 

Format  5E10.0 

RVHT  Reentry  altitude  (ft) 

(If  RVAT =0 . ,  the  default,  300000.  ft,  will  be  used) 

RVANG  Reentry  angle  at  altitude  RVHT  (deg) 

RVVEL  Reentry  velocity  at  altitude  RVHT  (ft/sec) 

RVMAS  Vehicle  mass  (slug) 

RVARA  Vehicle  reference  (base)  area  (ft2) 
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Card  10  Drag  Data  Array  Parameters 

(Enter  only  if  ITYPE=4,  RV,  and  NTRJ  <  0) 

Format  215 


Format 

NDG1 


Number  of  Mach  Number  values  input  on  card  data 
sets  as  Card  11  (Maximum  of  25) 

Number  of  angles  of  attack  input  on  card  data  sets 
as  Card  12  (Maximum  of  10) 


Card  11  Mach  Numbers  for  Drag  Coefficients 
-  (Enter  onlyTf  TTYPE-4.  RV."  ahd 


Format 


7E10.0 


Mach  number  values  in  increasing  magnitude 
(Enter  NOG!  values  to  maximum  of  25) 


Card  12  Angles  of  Attack  for  Drag  Coefficients 

(Enter  only  if  ITYpE=4,  RV,  and  NTRJ  <  0) 

Format  7E10.0 

ANGX  Angle  of  attack  values  in  increasing  magnitude.  The 

first  val ue  must  be  0. 

(Enter  NDG2  values  to  maximum  of  10) 

Card  1 3  Normal  Drag  Force  Coefficients 

(Enter  only  ITYPE=4,  RV,  and  NTRJ  <  0) 

Format  7E10.0 

CN  Normal  drag  force  coefficients.  Enter  NDG2  data  sets 

corresponding  to  the  angles  of  attack.  Each  data  sets 
contains  NDG1  values  corresponding  to  the  Mach  numbers 
(Total  number  of  values,  NDGlxNDG2,  maximum  in  number 
of  values,  25x10) 


Card  14  Axial  Drag  Force  Coefficients 


Format 


(Enter  only 
7E10.0 


~  and  NTRJ  <  0) 


CA  Axial  drag  force  coefficients.  Enter  NDG2  data  sets 

corresponding  to  the  angles  of  attack.  Each  data  sets 
contains  NDG1  values  corresponding  to  the  Mach  numbers. 
(Total  number  of  values,  NDGlxNDG2,  maximum  number  of 
values,  25x10) 

Card  15  Reentry  Vehicle  Geometry 

(Enter  only  if  ITYPE=4,  RV,  and  NTRJ  >  0) 

Format  7E10.0 

RVNR  Nose  radius  (in) 

RVBR  Base  radius  (in) 

RCVHA  Cone  half  angle  (deg) 

RVSR  Surface  roughness  (in) 

(If  RVSR=0.,  the  default,  0.02  in,  is  used) 
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RVCYHA  Cylinder  half  angle  (deg) 

RVFHA  Flare  half  angle  (deg) 

RVCL  Cone  length  along  axis  (in) 

RVFL  Flare  length  along  axis  (in) 

RVCYHA,  RVFHA,  RVCL,  RVCYL,  and  RVFL  are  not  presently  operational 
in  this  version  of  the  code,  but  blank  fields  at  least  are  required 
in  the  input  stream. 


Card  16  Booster  Trajectory  Data 

(Enter  only  if  ITYPE=2,  booster,  and  NTRJ  <  0) 

Format  5E10.0 


TIMX 

ALT 

RNG 

VEL 

ANG 

Enter  jNTRJl 


Time  after  launch  (sec) 

Altitude  (ft) 

Ground  range  (NMi) 

Velocity  (ft/sec) 

Flight  path  angle  to  the  local  horizon  (deg) 
card  data  sets  (Maximum  of  75) 


Card  17  Detonating  RV  Identifiers  and  Flags 

Format  A4,  IX,  A1  ,  4X,  15 

LTYPE2  Detonating  RV  (RVs)  identifiers 

LH0B  Height  of  burst  flag 

IMRV  MRV  flag 

0  -  Single  RV 

1  -  Circular  pattern  of  3  or  more  RVs 

2  -  Two  in-line  RVs 

3  -  Two  in-line  RVs  of  different  yield.  The  uprange 

burst  can  have  a  different  yield,  weapon  charac¬ 
teristics,  burst  altitude  and  be  at  a  later  time. 
(This  option  is  not  presently  implemented) 

Card  18  MRV  Parameters 

(Enter  only  if  IMRV#)) 

Format  2E10.0 

RADMAX  MRV  pattern  radius.  Half  the  ground  range  distance 
between  bursts  for  IMRV=2  or  3.  (NMi) 

RADTIM  Time  between  bursts  (sec)  (Enter  only  if  IKRV=3) 


Card  19  Detonating  RV  Burst  Altitude 
Format  El  0.0 

XH0B2  Burst  altitude  for  RV  (ft) 
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Card  20  Warhead  Data  for  Booster  Exclusion  Region  Calculations 
(Enter  only  if  ITVP£=2,  Booster) 

Format  2E10.0,  15,  E10.0,  15 

YLD  Warhead  yield  (kT) 

XNEUTX  Neutron  output  (n/kt) 

INEUT  Neutron  transmission  flag 

For  INEUT  >  0,  use  prestored  neutron  transmission  data  set  number  IXRAY 
and  if  XNEUTX=0.  use  corresponding  prestored  neutron  output.  For 
INEUT  <  0,  input  neutron  transmission  data  using  Cards  22  and  23  XNEUTX 
must  be  input. 

XFRAC  X-ray  energy  fraction 

IXRAY  X-ray  transmission  flag 

For  IXRAY  >  0  use  prestored  X-ray  transmission  data  set  number  IXRAY 
and  if  XFRAOO.  use  corresponding  prestored  x-ray  energy  fraction.  For 
IXRAY  <  0  input  x-ray  transmission  data  using  Cards  24  and  25  and  XFRAC 
must  be  input. 


Card  21  Warhead  Data  for  RV  exclusion  Region  Calculations 
(Enter  only  if  ITYP£=4,  RV) 

Format  2E10.0,  15,  E10.0,  15,  2E10.0 

YLD  Warhead  yield  (kt) 

XNEUTX  Neutron  output  (n/kt) 

INEUT  Neutron  transmission  flag 

For  INEUT  >  0  use  prestored  neutron  transmission  data  set  number  INEUT 
and  if  XNEUTXO.  use  corresponding  prestored  neutron  output.  For 
INEUT  <  0,  input  neutron  transmission  data  using  Cards  26  and  27  INEUTX 
must  be  input. 

XGAM  Gamma  ray  energy  fraction 

IGAM  Gamma  ray  transmission  flag 

GAMPL  Garrma  ray  pulse  width 

(If  GAMPL=Q.,  code  calculated  value  is  used) 

GAMDS  Gamma  ray  energy  fluence  to  dose  conversion  factor 
(rad(Si )-cm2/cal ) 

For  IGAM  >  0  use  prestored  gamma  ray  transmission  data  set  IGAM  and 
if  GAMDS=0.  use  corresponding  prestored  gamma  ray  energy  to  dose  con¬ 
version  factor.  For  IGAM  <  0,  input  gamma  ray  transmission  data  using 
Cards  26  and  27  GAMPL  must  be  input. 

Card  22  Neutron  Transmission  Data 
(Enter  only  if  INElIT  <  0) 

Format  7E1Q.0 

RH0RN  Air  mass  integral  values  for  neutron  transmission  data 
in  increasing  magnitude.  Enter  [INEUT]  values  (maximum 
of  60)  and  first  value  must  be  C.  (gm/cm2) 
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Card  23  Neutron  Transmission  Data 


XNETX  Normalized  neutron  transmission  values  corresponding  to 
the  air  mass  integral.  Enter  jlNEinj  values  (maximum 
of  60) 

Card  24  X-ray  Transmission  Data 

(Enter  only  if  iXRAY  <  0) 


Format 


7E10.0 


RH0RX  Air  mass  integral  values  for  x-ray  transmission  data 

in  increasing  magnitude.  Enter  JIXRAYJ  values  (maximum 
of  60)  and  first  value  must  be  0.  (gm/cm2) 

Card  25  X-ray  Transmission  Data 

(Enter  only  if  IXRAY- <  0) 

Format  7E10.0 

IPHI  Normalized  x-ray  transmission  values  corresponding  to  the 

air  mass  integral.  Enter  J IXRAY [  values  (maximum  of  60) 

Card  26  Gamma  Ray  Transmission  Data 


nter  only  v 
Format  7E10.0 

RH0RG  Air  mass  integral  values  for  gamma  ray  transmission  data 
in  increasing  magnitude.  Enter  jIGAMj  values  (maximum 
of  60)  and  first  value  must  be  0.  (gm/cm2) 


RH0RG 


Card  27  Gamma  Ray  Transmission  Data 
(Enter  only  if  IGAM  <  0) 

Format  7E10.0 

Normalized  gamma  ray  transmission  values  corresponding  to 
the  air  mass  integral.  Enter  /IGAM/  values  (maximum  of  60) 

For  IMRV=3  Card  19  and  Card  20  (booster,  ITYPE=2)  or  Card  21  (RV,ITYPE=4) 
are  input  and  Cards  22  through  27  are  input  as  necessary.  Three  addition 
conditions  for  this  pass  through  these  cards  exist.  (This  option  is  not 
currently  implemented). 

For  XNEUTXO.  and  INEUT=Q.,  the  neutron  transmission  data  and  output  are 
the  same' as  for  the  first  RV. 

For  XFRAC=0.  and  IXRAY=0,  the  x-ray  transmission  data  and  energy  fraction 
are  the  same  as  for  the  first  RV. 

For  GAMDS  and  IGAM=0,  the  gairnia  ray  transmission  data  and  energy  to  dose 
conversion  factor  are  the  same  as  for  the  first  RV. 


2.6  OUTPUT  DESCRIPTION 

The  output  of  the  NWEM  computer  code  has  been  structured  to  allow 
the  user  to  select  the  output  for  his  application  or  need  for  detail.  By 
allowing  this  flexibility,  the  NWEM  code  can  be  used  for  system  studies 
other  than  just  as  an  exclusion  region  contour  generator  for  AMM.  However, 
because  NWEM  was  developed  to  generate  exclusion  region  contour  data  for 
AMM,  the  output  will  obviously  reflect  this  aspect  of  the  NWEM  code. 

The  amount  of  output  the  user  chooses  is  controlled  by  the  input 
variable,  LEV0UT,  on  Input  Card  1.  This  choice  of  output  amount  or  level 
ranges  from  selecting  only  the  final  AMM  output  on  the  AMM  output  file 
(16)  to  selecting  detailed  input  echoing  and  detailed  exclusion  region  con¬ 
tour  data  by  environment.  In  particular,  there  are  four  levels  of  output. 
The  first  level,  LEV0UT  =  0,  is  only  that  data  for  the  final  AMM  exclusion 
region  contours  output  on  file  16  only.  The  second  level,  LEV0UT  =  1,  is 
the  additional  output  on  the  output  file  (6)  of  a  minimum  of  input  echoing, 
and  output  of  the  final  AMM  exclusion  region  contours.  The  third  level, 
LEV0UT  =  2,  of  output  increases  the  amount  of  data  output  to  the  output 
file  (6)  by  increasing  the  input  echoing  to  include  information  implied  by 
the  input  such  as  prestored  data  and  data  calculated  by  the  RV  trajectory 
routines  and  by  adding  the  AMM  exclusion  region  contour  data  by  environment. 
The  fourth  level,  LEV0UT  =  3,  adds  the  output  of  detailed  exclusion  region 
contour  data  for  several  environments.  For  all  levels  of  output,  error 

messages  and  other  significant  program  messages  are  always  output  on  the 

output  file  (6). 

As  stated,  the  objective  of  the  NWEM  computer  code  is  to  generate  AMM 
exclusion  region  contours  for  boosters  and  RVs  and  as  indicated,  this  infor¬ 
mation  is  always  output  on  output  file  16.  This  data  has  required  para¬ 
meter  orders  and  a  set  number  of  parameters  for  booth  boosters  and  RVs  to 
describe  the  exclusion  regions  treated  by  AMM.  These  regions  are  described 

in  detail  in  Volume  I  but  are  iterated  here  for  clarity.  The  exclusion 

region  contours  for  the  booster  are  a  static  circular  region,  a  dynamic 
segmented  annular  or  circular  "blast"  region  and  a  dynamic  rectangular 
region.  For  the  RV  exclusion  region  contours,  they  are  a  static  elliptical 
region  and  a  dynamic  rectangular  region.  A  listing  of  these  parameters. 
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their  order  and  their  computer  types  are  given  on  the  following  pages. 

This  list  also  corresponds  to  the  computer  variables  of  file  16  listed  at 
the  end  of  each  sample  problem.  For  additional  information  the  user  is 
referred  to  the  coimion  definitions  for  commons  / BST/  and  /RVC/. 

Because  the  output  to  the  normal  output  file  (6)  is  self  explanatory 
and  because  the  sample  problems  have  been  additionally  annotated,  the  user 
is  referred  to  the  sample  problems  for  explanation  of  the  output  to  this 
file. 
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AMM  Exclusion  Region  Contour  Parameters  and  Parameter 
Order  for  Booster  Flyout  Data 


Parame ter  Type 

Identification  of  booster  type  hollerith 

Identification  of  detonating  RV  type  hollerith 

Height  of  burst  option  for  RV  hollerith 

Radius  of  static  region  (NMi)  real 

Time  duration  of  static  region  (sec)  real 

First  time  for  blast  exclusion  region  (sec)  real 

Outer  radius  of  blast  exclusion  region  real 

at  first  time  (NMi) 

Inner  radius  of  blast  exclusion  region  real 

at  first  time  (NMi ) 

Angle  for  near  sector  of  blast  exclusion  region  real 

at  first  time  (deg) 

Angle  for  near  sector  of  blast  exclusion  region  real 

at  first  time  (deg) 

•  i 

•  • 

•  • 

Seventh  time  for  blast  exclusion  region  (sec)  real 

Outer  radius  of  blast  exclusion  region  real 

at  seventh  time  (NMi) 

Inner  radius  of  blast  exclusion  region  real 

at  seventh  time  (NMi) 

Angle  for  near  sector  of  blast  exclusion  region  real 

at  seventh  time  (deg) 

Angle  for  near  sector  of  blast  exclusion  region  real 

at  seventh  time  (deg) 

First  time  of  dynamic  region  (min)  real 

Downrange  dimension  of  dynamic  region  real 

at  fi rst  time  (NMi ) 

Cross  range  dimension  of  dynamic  region  real 

at  first  time  (NMi) 

Distance  from  burst  to  center  of  the  dynamic  region  real 

at  fi rst  time  (NMi ) 

•  • 

•  • 

•  • 

Tenth  time  of  dynamic  region  (min)  real 

Downrange  dimension  of  dynamic  region  real 

at  tenth  time  (NMi) 

Cross  range  dimension  of  dynamic  region  real 

at  tenth  time  (NMi) 

Distance  from  burst  to  center  of  the  dynamic  region  real 

at  tenth  time  (NMi ) 


AMM  Exclusion  Region  Contour  Parameters  and  Parameter 
Order  for  RV  Fratricide  Data 


i 


Parameter  Type 

Identification  of  incoming  RV  type  hollerith 

Identification  of  detonating  RV  type  hollerith 

Height  of  burst  option  for  detonating  RV  hollerith 

Reentry  angle  option  of  incoming  RV  (deg)  real 

Downrange  dimension  of  static  region  (NMi)  real 

Cross  range  dimension  of  static  region  (NMi)  real 

Uprange  dimension  of  static  region  (NMi)  real 

Time  duration  of  static  region  (sec)  real 

First  time  of  dynamic  reqion  (min)  real 

Half-length  dimension  of  dynamic  region  real 

at  first  time  (NMi ) 

Half-width  dimension  of  dynamic  region  real 

at  first  time  (NMi ) 

Distance  from  burst  to  center  of  dynamic  region  real 

at  first  time  (NMi) 

•  • 

•  • 

•  • 

Tenth  time  of  dynamic  region  (min)  real 

Half-length  dimension  of  dynamic  region  real 

at  tenth  time  (NMi ) 

Half-width  dimension  of  dynamic  region  real 

at  tenth  time  (NMi ) 

Distance  from  burst  to  center  of  dynamic  region  real 

at  tenth  time  (NMi ) 


1 


l 
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2.7  SAMPLE  PROBLEMS 


Two  sample  problems  demonstrating  the  capabilities  of  the  NWE M 
computer  code  are  presented  in  this  section.  These  problems  are  gener¬ 
ally  representative  of  boosters  and  RVs  and  do  not  represent  any  partic¬ 
ular  system 

The  first  sample  problem  is  a  booster  flyout  exclusion  problem. 

Input  Card  1  indicates  that  this  problem  is  a  booster  identified  as  BST1  , 
prestored  trajectory  data  set  1  is  used,  all  output  levels  are  printed 
and  the  problem  title  is  BOOSTER  ALL  LOW  ALTITUDE  ENVIRONMENTS.  Input 
Cards  2  through  8  indicate  a  1.0  psi  free  field  overpressure  is  the  only 
blast  criteria,  the  25  min.  cloud  cut-off  time  is  used  for  0.3  kj/cm: 
intercepted  kinetic  energy  dust  criteria,  cloud  avoidance  for  pebble/ejecta 
is  in  effect,  a  1.0  cal/cm2  thermal  fluence,  thermal  criteria,  is  used  with 
the  1.0  default  thermal  transmission  factor,  l.OxlO12  n/ cm:  neutron  fluence 
is  the  neutron  criteria,  a  0.1  cal/cm2  x-ray  fluence  is  the  x-ray  criteria 
and  EMP  is  not  considered.  Input  card  17  indicates  that  the  detonating 
RV  is  identified  as  RV  A  with  a  height  of  burst  option  of  A, and  it  is  not 
an  MRV  problem.  Card  19  indicates  a  burst  altitude  of  100.  ft.  and  Card  20 
indicates  that  the  yield  is  1000.  kt,  that  the  prestored  neutron  transmis¬ 
sion  data  set  1  and  its  respective  neutron  output  are  used,  and  that  the 
prestored  x-ray  transmission  data  set  2  and  its  respective  x-ray  energy 
fraction  are  used. 

The  second  sample  problem  is  an  RV  fratricide  exclusion  problem. 

Input  Card  1  indicates  that  this  is  an  RV  identified  as  RV  B,  that  the  drag 
data  are  calculated  as  well  as  the  trajectory,  all  output  levels  are 
printed  and  the  problem  title  is  ALL  REENTRY  ENVIRONMENTS.  Input  Cards  2 
through  4  indicate  that  a  250.  G  total  acceleration  is  the  only  blast  cri¬ 
teria,  that  the  25  min.  default  cloud  cut-off  time  is  used  for  the 
1.0  kJ/cm2  intercepted  kinetic  energy  dust  criteria  and  that  cloud  avoid¬ 
ance  for  pebble/ejecta  is  in  effect.  Card  6  indicates  a  1 . 0x1 0 1 2  n/cm2 
neutron  fluence  is  the  neutron  criteria  and  Card  7  indicates  a  peak  gamma 
ray  dose  rate  of  1 .0x1 0 7  rad(Si)/sec  is  the  gairnia  ray  criteria.  Card  9 
indicates  that  the  initial  reentry  altitude  is  300000.  ft,  initial  reentry 
angle  is  24°,  initial  reentry  velocity  is  24000.  ft/sec,  the  RV  mass  is 
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22.9  slugs  and  the  RV  reference  area  is  5.15  ft2.  Card  15  indicates 
that  the  RV  nose  radius  is  3.84  in.,  the  RV  base  radius  is  15.37  in., 
the  RV  cone  half  angle  is  7.5°,  the  0.02  in.  default  surface  roughness 
is  used  and  the  other  variables  ignored  but  necessary  as  explained  by  tne 
input  guide.  Input  Card  17  indicates  that  the  detonating  RV  is  identi¬ 
fied  as  RV  A  with  a  height  of  burst  option  of  S  and  it  is  not  an  MRV  pro¬ 
blem.  Card  19  indicates  a  burst  altitude  of  sea  level  and  Card  20 
indicates  that  the  yield  is  1000.  kt,  that  prestored  neutron  transmission 
set  1  and  its  respective  neutron  output  are  used  and  that  prestored  gamma 
ray  transmission  data  set  1  with  its  respective  gamma  ray  energy  to  dose 
conversion  factor  are  used  with  the  gamma  ray  energy  fraction  ana  pulse 
width  calculated  by  NWEM. 

From  the  expanded  output  for  these  cases,  it  is  verified  that  the 
inputs  to  the  problems  have  been  entered  and  stored  appropriately .  In 
addition,  the  selected  prestored  trajectory  data,  neutron  transmission 
data  and  x-ray  transmission  data  are  listed  and  can  be  verified  for  the 
booster  problem  (sample  problem  1  ),  and  the  calculated  drag  coeff *' c i cr f s , 
calculated  RV  trajectory  data,  prestored  neutron  transmission  data  and 
gamma  ray  transmission  data  are  listed  and  can  be  verified  for  the  RV 
problem  (sample  problem  2).  The  next  output  set  is  the  mechanical  environ 
men t  data  consisting  of  the  detailed  blast  exclusion  region  data  by  time 
increments,  followed  by  the  reduced  (AMM  region  type)  blast  exclusion  re¬ 
gion  data,  the  reduced  ejecta  exclusion  region  data  and  the  combined 
reduced  pebble/dust  exclusion  region  data.  Following  the  mechanical 
environment  output  is  the  output  of  the  radiation  environments  data,  con¬ 
sisting  of  detailed  radiation  environment  data  followed  by  the  reduced 
radiation  exclusion  region  data.  Following  the  radiation  environment 
data  is  the  printed  output  of  the  final  AMM  exclusion  region  data  wnich 
agrees  with  the  data  output  to  the  AMM  exclusion  region  output  file  (16). 

To  facilitate  user  understanding  of  the  output,  the  output  listings 
of  the  sample  problems  have  been  annotated  to  reflect  their  level  and 
def ini tion . 
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STORED  TRAJECTORY  DATA  USED 
SET  NUMDER  l 
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3.0  COMPUTER  CODE  MAINTENANCE  AND  UPDATE 


The  NWEM  computer  code  was  developed  with  a  philosophy  to  make 
maintenance  and  update  easy.  With  this  philosophy  in  mind,  the  code  was 
structured  in  a  modular  form  for  ease  of  changing  not  only  individual  sub¬ 
routines  or  functions  but  complete  environmental  model  and  exclusion  region 
contour  generating  routines,  of  extending  the  allowed  exclusion  region 
shapes  and  time  increments  and  of  expanding  the  prestored  data  base.  By 
having  the  flexibility,  the  NWEM  computer  code  can  virtually  be  modified 
to  perform  a  variety  of  studies  depending  on  the  complexity  of  the  en¬ 
vironmental  models  and  the  detail  of  the  exclusion  region  data  used.  The 
objective  of  this  section  is  not  to  define  studies  to  be  performed  by  NWEM 
but  to  give  the  user  information  necessary  to  make  modifications  to  the 
NWEM  computer  code  if  needed.  To  this  end,  this  section  discusses  data 
base  extension,  environmental  model  modification  and  extensions  of  the 
exclusion  regions. 

3.1  DATA  BASE  EXTENSION 

As  previously  stated  in  the  block  data  description  section,  it  was 
decided  to  prestore  input  quantities  in  block  data  that  are  needed  for 
the  radiation  transmission  data  and  booster  trajectory  data  to  facilitate 
code  check-out.  To  a  limit,  this  technique  of  prestored  data  can  be  ex¬ 
tended  to  include  increasingly  more  sets  of  transmission  data  and  tra¬ 
jectories  . 

If  a  user  desires  to  make  these  changes,  he  will  be  required  to  expand 
the  appropriate  dimension  statement  limits  designating  the  number  of 
sets  of  data  stored  and,  if  necessary,  the  dimension  statement  limits  de¬ 
signating  the  number  of  values  stored  per  data  set.  In  addition,  there 
are  variables  specifying  storage  limitations  which  may  require  change. 

As  the  code  is  presently  set  up,  it  is  relatively  easy  to  add  sets  of 
data  with  the  number  of  values  per  set  equal  to  or  less  than  the  present 
storage  limitations.  This  type  of  addition  of  data  is  recommended.  To 
change  the  numbers  of  data  values  stored  per  set  requires  updating  all 
common  blocks  associated  with  that  data  throughout  the  code,  but  will  not 
require  any  other  re programming  since  those  parts  of  the  code  using  this 
data  are  programmed  to  handle  a  variable  number  of  data  values. 


To  make  use  of  this  technique  for  booster  trajectories  for  adding 
data  sets  of  less  than  or  equal  values,  it  is  necessary  to  increase  the 
index  of  the  /TRAJX/  comnon  variables  in  block  data  and  subroutine  INPREP 
from  its  present  value  of  1  to  the  appropriate  number  and  enter,  in  the 
appropriate  order,  in  block  data  the  trajectory  data  as  a  function  of  time 
after  launch,  the  number  of  data  values  per  set  and  the  number  of  sets. 

If  there  is  a  change  required  of  the  first  index,  then  the  additional  change 
of  the  value  of  the  variable  NTRJM  (conmon  /TRNSL/ )  and  common  blocks  /TRAJ/ 
and  /HTRAJ/  must  be  changed  appropriately  in  each  subroutine  in  which  they 
appear. 

For  the  radiation  transmission  data  extension,  a  similar  procedure  is 
used.  Depending  on  the  environment  updated,  the  index  of  those  variables  in 
common  / S0R/  corresponding  to  that  environment  are  increased  and  values 
entered.  For  increases  in  the  number  of  values  per  set,  it  is  necessary 
to  change  the  value  of  the  variables  INETM  and  IXRM  appropriately,  and 
the  indexes  of  common  block  /0UT/  must  be  changed  appropriately  for  each 
subroutine  in  which  it  appears. 

As  mentioned  earlier  this  technique  of  data  extension  will  reach  a 
limit.  This  limit  is  the  restriction  computer  code  size  places  on  the 
user.  Since  these  dimensioned  variables  require  storage,  eventually  the 
core  size  will  be  exceeded  and  another  method  of  access  of  the  data  will 
be  necessary.  A  logical  choice  is  one  or  more  data  files.  As  the  NWEM 
computer  code  is  currently  written,  minor  modifications  are  required  to 
implement  this  method.  The  only  changes  would  be  to  replace  the  current 
block  data  logic  with  file  reading  logic  in  the  INPREP  subroutine. 

For  any  data  extensions  undertaken,  however,  the  user  is  cautioned 
to  review  in  detail  the  definitions  of  the  common  block  concerned  and 
be  very  familiar  with  the  INPREP  subroutine. 

3.2  ENVIRONMENT  MODEL  REPLACEMENT 

3ecause  the  NWEM  computer  code  has  been  structured  as  a  compilation 
of  independent  environment  models  and  exclusion  region  contour  generating 
routines,  replacements  or  exchange  of  environment  models  is  simplified. 
Additionally,  one  of  the  objectives  of  presenting  the  subroutine  descrip¬ 
tions  and  common  block  descriptions  in  the  formats  used  is  to  enable 


the  user  to  completely  understand  the  input  and  output  variables  of  each 
routine  and  to  trace  the  flow  of  the  calculation  down  through  the  other 
routines.  With  this  information,  the  user  can  determine  routines  and 
variables  affected  and  correctly  make  the  necessary  changes  to  supply  the 
rest  of  NWEM  with  the  appropriate  data  from  his  replacement  routines. 

Obviously,  it  would  be  misrepresentation  to  imply  that  environment 
model  replacement  can  be  undertaken  by  anyone  or  without  considerable  ex¬ 
perience  with  the  code.  However,  once  this  experience  has  been  gained, 
the  user  should  be  able  to  effect  the  changes. 

As  an  example,  assume  that  the  RV  trajectory  is  to  be  replaced. 

(While  this  example  is  not  an  environment  model,  it  will  serve  as  a  straight 
forward  example).  The  input  to  the  present  routines  is  the  RV  geometry 
parameters  (see  subroutine  PTMASS)  and  possibly  the  drag  force  coeffi¬ 
cients,  coirmon  /DRAG/,  and  the  output  is  the  RV  trajectory,  cotimon  /TRAJ/, 
and  possibly  again  the  drag  force  coefficients,  cotimon  /DRAG/.  To 
effect  the  change,  the  new  routines  only  need  to  be  able  to  use  this  input 
data  and  output  the  required  data  to  NWEM.  Those  subroutines  and  func¬ 
tions  (in  NWEM)  compatible  with  those  in  the  replacement  routines  can  be 
used  or  for  minimum  interactions  no  other  connections  between  NWEM  and 
the  replacement  routines  are  needed. 

3.3  EXTENSION  OF  ALLOWED  EXCLUSION  REGION  CONTOURS 

This  area  of  code  update  is  the  most  complex  of  the  three  considered. 
Because  there  exist  restrictions  imposed  by  the  limitations  and  restric¬ 
tions  placed  on  the  exclusion  region  contours  by  AMM,  routines  had  to  be 
developed  which  would  combine  the  individual  environmental  exclusion  rou¬ 
tines  and  construct  the  AMM  exclusion  region  contours  (see  Volume  1). 

To  extend  the  AMM  exclusion  region  contours,  the  user  must  be  very  famil¬ 
iar  with  the  function  of  these  special  routines.  However,  the  individual 
environment  models  and  exclusion  region  contour  generating  routines  pro¬ 
vide  detailed  exclusion  region  contour  information  which  can  be  used 
as  they  are  presently  as  a  bases  for  the  AMM  extended  exclusion  region 
contours.  The  combining  routines  modified  to  provide  the  desired  data. 

Since  any  realistic  example  of  extensions  of  this  type  are  quite  involved 


and  require  knowledge  of  the  code,  none  will  be  provided  in  the  hope  that 
before  the  user  attempts  it  he  will  be  familiar  with  the  code  and  will  be 
able  easily  to  effect  a  change  of  this  type. 

One  type  of  extension,  however,  that  only  requires  an  understanding 
of  the  code  is  to  use  the  detailed  exclusion  region  data  and  manipulate 
this  data  directly.  Essentially  all  special  routines  are  eliminated  and 
replaced  with  those  of  the  user's  design.  While  this  method  is  easier 
from  a  code  understanding  position,  it  may  require  considerable  code  de¬ 
velopment  effort  depending  on  the  required  type  of  output  data. 
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